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ZOOLOGICAL RESEARCH 


2018 New Year Address of Zoological Research 


At the beginning of a wonderful new year, we look back at 
Zoological Research (ZH) in 2017. We are very grateful to 
all our readers and authors for your dedication and continued 
support of ZH. Your ideas, input, and enthusiasm have been of 
immense value in helping us to improve the journal. Here, we 
would like to share a few memorable events and people of the 
past year. 

Firstly, it is with deep regret and sorrow that we note 
the passing of Colin Peter Groves, Professor of Biological 
Anthropology from the Australian National University in 
Canberra, Australia, on 30 November 2017. He was a 
long-standing and loyal friend of ZH, and joined our editorial 
board in 2014. He fulfilled his duties with dedication and 
diligence, and shared his outstanding expertise with opinions 
and constructive suggestions on many submissions to the 
journal. Despite his poor health, he continued to contribute 
to the journal (Groves, 2016) and was a respected member 
of our editorial board. We greatly cherish his contributions 
and generosity. 

We have also been inspired in the past year by the continued 
improvement in academic quality, total citations, and general 
influence of each of our publications. ZR achieved an 
estimated impact factor of 0.73 (based on the citation data from 
the Web of Science) and CiteScore of 0.84 (dated 11 January 
2018, Scopus). We have also been fortunate with the addition 
of eleven outstanding experts who recently joined ZH as 
editorial members, including Yu-Hai Bi (Institute of Microbiology, 
Chinese Academy of Sciences (CAS), China), Peng-Fei Fan 
(Sun Yat-Sen University, China), Patrick Giraudoux (University 
of Franche-Comté, France), Cyril C. Grueter (University of 
Western Australia, Australia), Wei-Zhi Ji (Kunming Institute 
of Zoology, CAS, China), Shu-Qiang Li (Institute of Zoology, 
CAS, China), Wen-Jun Liu (Institute of Microbiology, CAS, 
China), Julian Kerbis Peterhans (Roosevelt University, USA), 
Xiang-Guo Qiu (University of Manitoba, Canada), Rui-Chang 
Quan (Xishuangbanna Tropical Botanical Garden, CAS, China), 
and Christian Roos (Leibniz-Institute for Primate Research, 
Germany). These remarkable specialists will bring new ideas 
and expertise to help reshape the journal. We are confident 
that our excellent editorial board will further facilitate the pivotal 
and active role of ZR in the field of science and publishing. 

In 2017, the "Project for Enhancing International Impact 
of China STM Journals" (PIIJ) (Class B) (2016-2018), the 
most extensive and influential journal fostering project in China, 
continued its support of ZR. With great effort from all editors 
and staff from the editorial office, ZR has been ranked among 
the top 300 "Outstanding S&T Journals of China". Moreover, 
with the successive release of special issues focusing on 
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academic topics of current interest, ZR has evolved into a 
vibrant journal with appreciable readability. 

In February 2017, ZH successfully hosted the "2017 
Frontiers in Zoology Symposium" with the theme of "Animal 
Genomics and Ecological Protection". We hope to replicate 
this success with the "2018 Frontiers in Zoology Symposium", 
which will focus on "Comprehensive Scientific Investigation of 
Animals on the Tibetan Plateau and East Asia" (please find 
details at www.zoores.ac.cn). As such, we hope to see you 
in Kunming this coming spring. ZH also co-hosted the 2017 
Annual Conference of the Chinese Herpetology Society. Such 
events are wonderful ways in which to meet up the old friends 
and colleagues as well as establish new contacts and share 
achievements and opinions. 

To facilitate scientific communication and promote subject 
development, ZH is not only increasing its academic 
value, but also consistently reinforcing the ethical integrity 
accompanying academic publications. In 2016, ZR published 
an editorial (Liu, 2016), letter to the editor (Joob & 
Wiwanitkit, 2016), and statement (http://www.zoores.ac.cn/ 
EN/abstract/abstract3772.shtml) regarding proper authorship. 
Publishing is of importance in almost every stage of a 
researcher's career (Editorial Office of Zoological Research, 
2016). Nowadays, under the background of big science, 
the advancement and elucidation of scientific questions often 
require collaboration among independent research groups, 
even those from different fields, to combine their expertise 
or specialties. | Academia benefits from such extensive 
collaborations. As a result, we have observed that the 
proportion of co-first or co-corresponding author articles has 
significantly increased during the last few decades in many 
scientific journals. Akhabue & Lautenbach (2010) studied 
original research articles with equally contributing authors 
among five top life science journals, i.e., New England Journal 
of Medicine, Journal of American Medical Association, Annals 
of Internal Medicine, Lancet, and British Medical Journal. In 
2000, the proportions of co-first or co-corresponding author 
papers in the five journals were «196, 0%, 096, «196, and 0%, 
whereas in 2009, the proportions had significantly increased 
to 4.4%, 2.396, 1.6%, 2.7%, and «196, respectively. Li et al. 
(2013) studied the numbers of equally contributing articles in 
three major anesthesia journals (Anaesthesia, British Journal 
of Anaesthesia, and Anesthesia & Analgesia) over a 10-year 
period, and found that such papers increased significantly from 
2002 to 2011 (0.9% to 8.8%, 0% to 8.8%, and 0.3% to 3.4%, 
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respectively). The Journal of Genetics & Genomics published 
124 and 85 original articles in 2007 and 2016, respectively, 
with co-first or co-corresponding author papers accounting 
for 1696 (20/124) and 41.196 (35/85), respectively. In Cell 
Stem Cell, the percentage of co-first and co-corresponding 
author articles increased from 20.5% (18/88) in 2007 to 37.6% 
(77/205) in 2016. This trend reinforces collaborations and data 
sharing. For ZH, in 2007, approximately 996 (9/108) of the 
original research articles published were by two corresponding 
or co-first authors; by 2016, this increased to 3396 (12/36). 

For most academic journals there are no universal rules 
regarding shared authorship, with different journals usually 
having their own specific regulations. For example, the American 
Journal of Human Genetics states that "unrestricted joint 
authorship is allowed" and “a maximum of two corresponding 
authors is allowed" (http://www.cell.com/ajhg/authors). Neuron 
merely restricts the number of lead contacts (i.e., the 
corresponding author responsible for communicating with 
the journal and helping facilitate other issues regarding 
submission) to only one, but not other authors (http: 
/www.cell.com/neuron/authors#policies). The Journal of 
Neurosciences does not restrict co-authors numbers, but 
states “honorary authorship” to be a “misrepresentation” (http:// 
www.sfn.org/member-center/professional-conduct/guidelines- 
for-responsible-conduct-regarding-scientific-communication# 
menulevel1). Scientific Reports allows up to six co-authors, 
with corresponding authors usually limited to three. 

ZR has an explicit attitude and firm stance to proper 
authorship and resolutely complies with the criteria of 
the International Committee of Medical Journal Editors 
(International Committee of Medical Journal Editors, 2013). 
Accordingly, ZR would consider guest/ghost authors unethical. 
ZR appreciates and encourages positive and comprehensive 
collaboration. Combining findings from different researchers 
can yield a more complete story, and emphasizing equal 
contribution of involved parties is fair and encourages team 
work and collaborations. However, it is commonly accepted, 
at least in the biomedical field, that publication credit can 
never be enjoyed equally among all listed authors in the 
byline, with the first and last (corresponding) authors usually 
receiving most recognition, unless shared authors are clearly 
listed in alphabetical order. The addition of equally contributing 
authors can dilute the credit allocated to each author. These 
actions obviously deviate from the original intention for properly 
sharing credit. 

Although ZR encourages every author of a paper to 
read the authorship guidelines thoroughly and carefully 
before submission, we have not previously insisted upon 
an author-contribution statement. However, we believe it is 
important to increase transparency and enhance authorship 
integrity. Thus, for submissions to ZH after 1 January 2018, 
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an author-contribution statement will be compulsory. The 
statement is not required in uniform format, but the contribution 
of every person in the author list must be specifically described. 
Moreover, for collaborative studies, ZR requires that every 
person involved in the paper agrees with the author list and 
sharing of the credit. 

Whether these measures will make a substantial difference 
in eliminating unethical authorship remains to be seen. 
However, we believe such ambitions are worth attempting 
simply because transparency and fairness are always worthy 
of pursuit. With your support and positive contributions, ZH will 
remain a respectable platform for scientific communication. 

Again, thank you all. May your new year hold all the 
excitement, promise, and good wishes we have shared in 2017. 


Sincerely, 


Yong-Gang Yao, Editor-in-Chief 
Kunming Institute of Zoology, 
Chinese Academy of Sciences, Kunming Yunnan 650223, China 


J ae: om 


rag 


Xue-Long Jiang, Executive Editor-in-Chief 
Kunming Institute of Zoology, 
Chinese Academy of Sciences, Kunming Yunnan 650223, China 
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ABSTRACT 


Due to their inability to generate a complete immune 
response, mice knockout for type | interferon (IFN) 
receptors (/fnar^) are more susceptible to viral 
infections, and are thus commonly used for 
pathogenesis studies. This mouse model has been 
used to study many diseases caused by highly 
pathogenic viruses from many families, including the 
Flaviviridae, Filoviridae, Arenaviridae, Bunyaviridae, 
Henipaviridae, and Togaviridae. In this review, we 
summarize the findings from these animal studies, 
and discuss the pros and cons of using this model 
versus other known methods for studying pathogenesis 
in animals. 


Keywords: lfnar; Mice; Animal model; Flavivirus; 
Filovirus; Arenavirus, Bunyavirus,  Henipavirus; 
Togavirus 


INTRODUCTION 


Outbreaks of infectious diseases amongst the human 
population have been documented for thousands of years. The 
earliest on record was the Plague of Athens between 429-426 
B.C.. Caused by an unknown pathogen, the outbreak killed 
over 75 000 people (Littman, 2009). Epidemics that have 
occurred since are too numerous to list thoroughly in this review, 
but include multiple instances of plague (Yersinia pestis) that 
devastated populations of Europe, Asia and North Africa with 
hundreds of millions of deaths. The most infamous epidemic of 
plague was the “Black Death" during 1346—1350 (Haensch et 
al., 2010), in which an estimated 3096—6096 of the population 
was wiped out. Viral outbreaks including those caused by 
smallpox, measles and viral hemorrhagic fevers in various 
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locations worldwide have impacted tens of millions (CDC 2017; 
Moss & Griffin, 2012; Théves et al., 2016). In the 20" century, 
three major influenza pandemics (H1N1 during 1918, H2N2 
during 1957-1958 and H3N2 during 1968-1969) have killed 
over 75 million people combined (Johnson & Mueller, 2002; 
Kilbourne, 2006). As we enter the 21*' century, the death toll 
from outbreaks of infectious diseases has decreased 
dramatically, and the highest numbers of fatalities were from the 
2014-2016 Ebola virus outbreak in West Africa (over 11 000 
deaths) (WHO, 2016a) and the 2009 H1N1 influenza pandemic 
(over 18 000 deaths) (WHO, 2010). Considerable advances 
and deployment of prophylactics, therapeutics, rapid point-of- 
care diagnostics and surveillance have limited the negative 
impacts from outbreaks in many parts of the world and saved 
many lives that would otherwise have been lost. However, 
outbreaks of re-emerging infectious diseases have been 
occurring with ever increasing frequency in recent years, and 
there is still much to do in the war against infectious diseases. 
The use of animals to study pathogenesis, as well as test 
potential vaccines and drugs, have played a big role in 
accelerating the most promising compounds through the pre- 
clinical process before testing in clinical trials. Non-human 
primates (NHPs), the closest relative species to humans, are 
considered the gold standard animal model for many infectious 
diseases (Safronetz et al., 2013) because these animals 
recapitulate multiple aspects of human disease, and thus any 
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experimental results are expected to have high translatability 
and applicability to humans. However, NHPs are costly to 
acquire, difficult to handle, and require specialized facilities to 
house and provide husbandry (Coleman, 2011), and are 
unaffordable for many research laboratories. To address this, 
smaller animals have been used for many preliminary studies 
and screens of candidate vaccines and drugs. For instance, the 
domestic ferret (Mustela putorius furo) was used to study many 
member viruses belonging to the order Mononegavirales 
(Enkirch & von Messling, 2015). 

Mice are an ideal species for studying human infectious 
diseases. The immune systems of mice and humans are often 
sufficiently similar that they can be infected with the same 
pathogens (Buer & Balling, 2003). For example, immunocompetent 
wild-type mice are susceptible to infections with a number of 
influenza virus subtypes (Belser et al., 2010; Driskell et al., 
2010; Gubareva et al., 1998; Xu et al., 2013), severe acute 
respiratory syndrome coronavirus (SARS-CoV) (Channappanavar 
et al., 2016) and Rift Valley fever virus (RVFV) (Smith et al., 
2010), and outbreaks with these pathogens can be rapidly and 
easily studied. Unfortunately, wild-type mice are not susceptible 
to many other pathogens with outbreak potential, and thus 
alternative strategies are needed. Mice lacking the type | 
interferon (IFN) receptor (/fnar^) were generated in 1994 
(Muller et al., 1994). While these transgenic mice do not show 
any overt abnormalities by six months of age and are fertile, the 
animals are entirely unresponsive to the effects of type | IFNs. 
Ablated immune responses in Ifnar^ mice were observed after 
challenge with Vesicular stomatitis virus, Semliki Forest virus, 
vaccinia virus, or lymphocytic choriomeningitis virus, and the 
knockout animals showed enhanced susceptibility resulting in 
either higher viral organ titers, or death at lower doses compared 
to wild-type mice (Muller et al., 1994). 

Since the type I IFN response plays such an important role in 
innate and adaptive immunity against viral infections (McNab et 
al, 2015), the /fnar^ mice, which are available in many 
backgrounds, have since been used to study many highly 
pathogenic viruses. In this review, we summarize the results of 
using /fnar^ mice to study selected pathogens from the 
Flaviviridae, Filoviridae, Arenaviridae, Bunyaviridae, Henipaviridae, 
and Togaviridae families, focusing on member viruses that 
either have, or may have the potential to cause large scale 
outbreaks in the future. 


FLAVIVIRIDAE 


The family Flaviviridae contains many member viruses which 
are highly pathogenic to humans and/or have high outbreak 
potential. West Nile fever, Dengue fever, Yellow fever, Japanese 
encephalitis, and Zika fever are all mosquito-borne diseases, 
caused by West Nile virus (WNV), Dengue virus (DENV), 
Yellow fever virus (YFV), Japanese encephalitis virus (JEV) and 
Zika virus (ZIKV), respectively. 


West Nile virus (WNV) 


Found in temperate and tropical regions, WNV is maintained in 
a mosquito-bird-mosquito cycle in nature, with humans as 
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incidental hosts. WNV was firstly identified in Uganda in 1937, 
but the majority of infections (~80%) caused only mild disease 
or were asymptomatic (WHO, 2011). In case of symptomatic 
disease, fever, headache, fatigue, muscle pain, nausea, 
vomiting and rash is observed. Less than 1% of cases are 
neuroinvasive, in which patients present high fever, stupor, 
disorientation, coma, tremors, convulsions, muscle weakness, 
and paralysis (WHO, 2011). Large outbreaks of WNV occurred 
sporadically throughout the decades, but during 1996 WNV re- 
emerged in Romania and caused 393 confirmed infections, in 
which 352 patients had manifestations in the central nervous 
system (Tsai et al., 1998). Subsequent epidemics of WNV with 
high rates of neuroinvasive disease was then noted in Morocco in 
1996, Tunisia in 1997, and large outbreaks in Italy and Israel in 
1998 (Hubálek & Halouzka, 1999). At present, WNV is endemic 
in Africa, Asia, Europe, Australia, and has spread into Canada 
and the United States (Chen et al., 2013). An outbreak of WNV in 
2012 in the United States claimed 286 lives (Murray et al., 2013). 

Wild-type 129Sv/Ev mice infected subcutaneous (SC) with 
10? plaque forming units (pfu) of WNV showed 62% mortality 
and died at a mean time to death of 11.9+1.9 days post 
infection (dpi), with no clinical signs observed until 8 dpi 
(Samuel & Diamond, 2005). In their study, Samuel et al. 
challenged 8—10 week old Ifnar^ mice (129Sv/Ev background) 
with 10°, 101 or 10? pfu or WNV strain 3000.0259 via footpad 
(SC) inoculation. Regardless of dose, all mice showed severe 
clinical symptoms by 3 dpi, including hunched posture, ruffled 
fur and reduced activities. Death occurred within 12-48 hours 
after the onset of symptoms, and the mean time to death was 
3.840.5 dpi for Ifnar~ mice in the 10? pfu group (Samuel & 
Diamond, 2005). Live infectious virus can be found in the 
muscle, heart, lung, kidney, liver, but not in the pancreas. 
Additionally, challenge of /fnar~ mice (C57BL/6 background) 
also showed 100% lethality and a mean time to death of 
3.4+0.5 dpi (Samuel & Diamond, 2005). 


Dengue virus (DENV) 

DENV is widespread in the temperate and tropical regions of 
the world, and each year approximately 50-150 million people 
are infected (Bhatt et al., 2013), with over 10 000 deaths 
(Stanaway et al., 2016). Symptoms of Dengue fever include a 
high fever, headache, vomiting, muscle and joint pains, and skin 
rash. Severe cases of disease is usually associated with 
secondary infection with heterologous types of DENV (Halstead, 
1988), and can develop into Dengue hemorrhagic fever (with 
hemorrhage, thrombocytopenia and blood plasma leakage), or 
into Dengue shock syndrome, both of which are potentially fatal 
(Kularatne, 2015). 

Infection of 129Sv/Ev mice with 109 pfu DENV-2 via the 
intravenous (IV) route resulted in 87% survival (26 out of 30 mice), 
and inoculation with 4.4x10^ pfu DENV-1 via IV resulted in 93% 
survival (40 out of 43 mice) (Shresta et al., 2004). In their study, 
Shresta et al. challenged 5-6 week old /fnar"- mice (129Sv/Ev 
background) with DENV-2 strain PLO46 (n=12) and DENV-1 
strain Mochizuki (n=16) at the same doses and inoculation 
routes. While no lethality was observed, sera and major organs 
harvested from infected mice at 3 and 7 dpi showed the 


presence of virus in all sera, liver, spleen, and lymph node 
samples, as well as some brain and spinal cord samples 
(Shresta et al., 2004). Interestingly, mice deficient for the type | 
and || IFN receptors (AG129) showed uniform death from 
DENV-2 infection with animals dying between 7—30 dpi, and 
DENV-1 as well with mice succumbing to disease between 7— 
14 dpi (Shresta et al., 2004). 


Yellow fever virus (YFV) 

YFV is endemic in tropical areas of Africa and South America 
(WHO, 2016b), when the virus was introduced via the slave 
trade during the 17" century. Many infections are symptomatic, 
but if clinical symptoms appear, they include fever, chills, 
appetite loss, nausea, muscle pains, and headaches. A small 
percentage (71596) of cases will go on to develop more severe 
disease including jaundice, dark urine, vomiting and abdominal 
pain. Hemorrhage from the mouth, nose, eyes or stomach may 
occur and 5096 of patients with these symptoms succumb to 
disease (WHO, 2016b). YFV was responsible for ~127 000 
severe infections and 45 000 deaths in 2013 (WHO, 2016b), 
with increased incidence over the past decades, and the risk of 
an outbreak in urban centers is a serious public health threat 
(Barrett & Higgs, 2007). 

Inoculation of wild-type 129 mice SC in each rear footpad 
with 10* pfu of YFV did not result in any weight loss or death 
(Meier et al., 2009). In their study, Meier et al. (2009) 
challenged 3—4 week old /fnar"- mice (129 background) with 
YFV strains Asibi or Angola73 under the same conditions. The 
mice were shown to be susceptible to the challenge, with death 
occurring between 7—9 dpi. Additionally, the mice developed 
viscerotropic disease with virus dissemination to the visceral 
organs, spleen and liver, in which severe damage of the organs 
can be observed with gross pathological examination and 
hematoxylin/ eosin staining. Elevated levels of MCP-1 and IL-6 
in these organs are suggestive of a cytokine storm (Meier et al., 
2009). 


Japanese encephalitis virus (JEV) 

Japanese encephalitis is an acute disease of the central 
nervous system in humans. Infected patients develop a febrile 
illness with headaches vomiting and diarrhea, as well as 
reduced levels of consciousness, seizures, and photophobia. 
Severe encephalitis occurs later in the disease course and is 
associated with a higher frequency of seizures, resulting in 
coma and death (Ghosh & Basu, 2009). Mental retardation may 
develop in the patient. JEV is endemic in large parts of Asia and 
the Pacific, and 30 000-50 000 infections (Solomon, 2006), 
including up to 15 000 deaths (Ghosh & Basu, 2009), are 
reported yearly. It is estimated that approximately 2596—309^6 of 
cases are fatal, but 5096 result in permanent neurological 
sequelae (Ghosh & Basu, 2009). 

Inoculation of wild-type 129 mice via the SC route with 10°, 
102, 10^ or 109 pfu of JEV strain JaOArS982 resulted in 
between 1096-4096 survival, but the difference was not 
statistically significant and deaths were not dose-dependent 
(Aoki et al., 2014). In their study, Aoki et al. repeated the SC 
infection in 5-6 week old /fnar~ mice (129 background) at the 


same doses. The mice were found to be very susceptible to the 
challenge, with uniform, dose-dependent death occurring at 64, 
80, and 96 hours after infection for the 109, 10* and 10? pfu 
groups, respectively. In the 10° group, 90% mortality rate was 
observed, and the animals died at 120 hours after infection. 
Live JEV could be detected in the spleens and brains of 
infected animals, with peak titers at 48 hours (Aoki et al., 2014). 


Zika virus (ZIKV) 
First isolated in 1947 from an infected monkey in Uganda and 
re-isolated from Aedes mosquitoes in the same area during 
1948 (Dick et al, 1952), ZIKV infections in humans have 
sporadically occurred in Africa and Asia, but in 2007 the virus 
continued spreading, causing outbreaks in small island 
countries located in the Pacific Ocean, such as Yap Island 
(Duffy et al., 2009), French Polynesia (Cao-Lormeau et al., 
2014) and Easter Island (Tognarelli et al., 2016). In early 2015, 
an epidemic of ZIKV infections, originating from Brazil, spread 
through most of North and South America and the Caribbean 
with tens of thousands of people over 80 countries infected 
(WHO, 2017), as well as thousands of imported cases from 
travelers returning to their home countries after visiting outbreak 
areas. The epidemic was declared over by the World Health 
Organization (WHO) on November 2016 (WHO, 2017), but 
many countries are still dealing with the long-term impact of 
ZIKV infections. Infections of ZIKV are typically asymptomatic, 
but if present they are mild in nature and includes fever, joint 
pain, maculopapular rash, and bloodshot eyes (Simpson, 1964). 
While no deaths have been reported from ZIKV infections, 
mother-to-child transmission during pregnancy may result in 
congenital Zika syndrome with abnormalities in the central 
nervous system (microcephaly, intellectual development, seizures 
and vision impairment) (Boeuf et al., 2016). ZIKV infections in 
adults is associated with Guillain-Barré syndrome (Frontera & 
da Silva, 2016). Distinct from other flavivirus infections, sexual 
transmission of ZIKV from male-to-male (Deckard et al., 2016), 
male-to-female (D'Ortenzio et al., 2016; Hills et al., 2016) and 
female-to-male (Davidson et al., 2016) have been documented. 
Infection of wild-type 129Sv/Ev mice SC with 109 pfu of ZIKV 
MP1751 did not result in any observable clinical symptoms or 
histological changes, despite the virus being detected at low 
levels in the blood, spleen and ovaries (Dowall et al., 2016). In 
their study, Dowall et al. challenged 5—6 week old Ifnar^ mice 
(129Sv/Ev background) under the same conditions as above, 
and showed that all animals succumbed to disease at 6 dpi with 
2096 body weight loss. High levels of virus could be detected 
by RT-qPCR at 3 and 7 dpi in the blood, spleen, brain, ovary 
and livers of these animals. Pathology studies show that 
inflammatory as well as degenerative changes could be seen in 
the brains of infected Ifnar^ mice (Dowall et al., 2016). In 
another study, Lazear et al. (2016) inoculated 5-6 week old 
Ifnar*~ mice (C57BL/6 background) with 10? pfu of ZIKV strain 
H/PF/2013 or MR766 via the SC route in the footpad. The 
results show that Ifnar^ mice all died within 8-10 dpi after 
challenge with H/PF/2013, and 80% death with MR766, with 
death between 9-13 dpi. Additionally, an SC challenge with 10? 
focus forming units (ffu) of ZIKV strain Dakar 41671, 41667 or 
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41519 in /fnar^- mice results in uniform death by 6 dpi (Lazear 
et al., 2016). In a third study, Rossi et al. (2016) inoculated 3-, 
5- and 11-week old /fnar^mice (C57BL/6 background) with 
1x10° pfu of ZIKV FSS13025 via the SC route. The results 
showed 10096 lethality in 3-week old animals with death 
occurring at 6—7 dpi, but only 5096 death in 5-week old animals 
and no deaths in 11-week old animals (Rossi et al., 2016), 
indicating that the disease caused by ZIKV infection in these 
animals is age-dependent. 


FILOVIRIDAE 


Ebola virus (EBOV), Sudan virus (SUDV), Reston virus 
(RESTV), Tai Forest virus (TAFV), Marburg virus (MARV) 
and Ravn virus (RAVV) 
The family Filoviridae consists of many member viruses, including 
EBOV, SUDV, RESTV, TAFV and MARV, among others. With the 
exception of RESTV, all filoviruses are pathogenic in humans and 
infected patients initially present with fever, sore throat, muscular 
pain, headaches, vomiting, and diarhea. As the infection 
develops, a rash is observed along with decreased organ 
function (especially liver and kidneys). Hemorrhage, shock and 
eventually multiple organ failure results in the death of the patient 
(Bradfute et al., 2012). Outbreaks of filovirus disease in humans 
are sporadic and unpredictable, and typically localized 
geographically to sub-Saharan Africa, but imported cases have 
occurred in the past to Europe and North America (CDC, 2014, 
2017b). The case fatality rate (CFR) of EBOV and MARV can 
reach up to 90%, whereas SUDV is ~50%. Only one case of 
TAFV has been recorded, in which the patient fell ill but survived 
infection (Formenty et al., 1999). The CFR of RAVV cannot be 
estimated accurately since the only large scale outbreak during 
1998-2000 in the Democratic Republic of the Congo (128 
deaths out of 154 cases) was due to the simultaneous co- 
circulation of RAVV and MARV (Bausch et al., 2003, 2006). 
Infection of wild-type adult immunocompetent 129 mice does 
not result in disease or death (Bray, 2001). Bray (2001) then 
inoculated 8—16 week old Ifnar^ mice (129 background) with 
1 000 pfu of EBOV, SUDV, RESTV, TAFV as well as MARV via 
intraperitoneal (IP) route. The results show that /fnar mice 
succumbed to infection with the Mayinga isolate of EBOV, with 
a mean time to death of 5.4 dpi, but resistant to the Kikwit 
isolate of EBOV. SUDV strain Boneface produced uniformly 
lethal infection with a mean time to death of 6.3 dpi, but RESTV 
and TAFV infections did not result in death of the Ifnar^ mice. 
Infection with RAVV and MARV produced 100% and 67% lethal 
infections, with a mean time to death of 6.0 and 8.5 dpi, 
respectively. Additionally, a SC challenge of the Mayinga isolate 
of EBOV to /fnar^- mice was shown to be fully lethal with a 
mean time to death of 7.3 dpi (Bray, 2001). In another study, 6— 
9 week old /fnar^- mice (129 background) were challenged with 
an aerosol dose of MARV between 10295? 50% tissue culture 
infective doses (TCIDso), EBOV between 10°? TCIDso of EBOV, 
or SUDV at 10*? TCIDa;. All animals challenged with MARV 
succumbed to disease at a mean time to death of 11—13 dpi, 
whereas EBOV-infected mice died at an average of 8 dpi. 
Clinical symptoms such as lethargy, weight loss and 
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piloerection were observed prior to death. Although symptoms 
such as anorexia were observed in SUDV-infected mice from 
7—11 dpi, all infected mice survived and returned to their pre- 
challenge weights by the conclusion of the experiment (Lever et 
al., 2012). 


ARENAVIRIDAE 


Lassa virus (LASV) 

Lassa fever is prevalent in the West African countries of 
Nigeria, Liberia, Sierra, Leone, Mali, Ghana, and Guinea, in which 
300 000—500 000 cases are reported yearly, including 5 000 
deaths (CFR ~1%) (Ogbu et al., 2007); however, the CFR from 
nosocomial outbreaks can reach as high as 65% (Fisher-Hoch 
et al., 1995). Carried by the multimammate rat (Mastomys 
natalensis), most infected patients are asymptomatic, but if 
illness occurs the initial presentation includes fever, weakness, 
headaches, vomiting, and muscle pains. In advanced disease, 
haemorrhaging, encephalopathy, shock and organ failure is 
observed (Schmitz et al., 2002). 

Wild-type 129S1SvlmJ mice are naturally resistant to 
infection with LASV (Yun et al., 2012). In their study, Rieger et 
al. infected 8-12 week old /fnar~ mice (129/Sv background) 
with 10? ffu of LASV strains Josiah, AV, BA366 and Nig04-10 
via the IV route. No deaths were observed with the mice, but 
the peak of viremia (10*5* ffu/mL of blood) was detected at 8 
dpi and still not fully cleared by 21 dpi. Weight loss of 
approximately 1596 by 8 dpi was observed, along with elevation 
of liver enzymes AST and ALT. Other findings include the 
presence of high levels (up to 10" ffu/g of tissue) of live LASV in 
the lung, kidney, heart, spleen, brain and liver infected animals at 
9-10 dpi (Rieger et al., 2013). Thus, the results support the 
establishment of productive LASV infection in Ifnar^ mice. 


BUNYAVIRIDAE 


Crimean-Congo hemorrhagic fever virus (CCHFV) 
Crimean-Congo hemorrhagic fever is caused by infections with 
CCHFV, which was first reported in the 1940s, but a study 
suggests that the virus may have been present since 1500- 
1100 B.C. (Carroll et al., 2010). Initial symptoms of CCHFV 
infection include fever (over 39.0 °C), muscle pains, fatigue, 
dizziness, vomiting, and diarrhea (Whitehouse, 2004). Advanced 
CCHFV infections are characterized by more severe symptoms 
including liver failure, petechiae as well as gastrointestinal and 
cerebral hemorrhage resulting in death (Whitehouse, 2004). 
The CFR can vary widely: it was reported to be 596 during an 
outbreak in Turkey (Kubar et al., 2011), but 6096 during another 
outbreak in the UAE (Schwarz et al, 1996). Spread by 
Hyalomma ticks, cases of CCHFV infections in humans has 
been reported in western Asia, Eastern Europe, the Middle East, 
as well as South Africa, although the geographical distribution 
of the Hyalomma vector is widespread and encompasses all of 
Africa, as well as European and Asian regions south of the N50” 
latitude (WHO, 2008). Approximately 50 cases are reported per 
year worldwide, but over 200 cases were reported during 2003— 
2004 (Messina et al., 2015). 


Infection of wild-type 129 Sv/Ew mice with CCHFV at high 
doses results in the establishment of an infection that is rapidly 
cleared from the kidney, brain, heart and blood within 3 dpi, and 
clearance from the liver and spleen by 11 dpi, but no clinical 
signs or mortality (Bereczky et al., 2010). In their study, 
Bereczky et al. challenged 7—10 week old /fnar^- mice (129 
Sv/Ew background) IP with 10'—109 ffu per animal of CCHFV 
strain IbAr 2000. Symptoms including laboured breathing were 
observed between 42-70 hours after infection and uniform 
death was observed at all doses. The highest viral loads in 
Ifnar^ mice were observed at 2 dpi in the spleen and liver (over 
10'? viral RNA copies/g of tissue), but could also be detected in 
blood, as well as other major organs including the kidney, brain 
and heart (Bereczky et al., 2010). In another study, Zivcec et al. 
infected 6—12 week old Ifnar^ mice (C57BL/6 background) with 
104^ TCIDso of CCHFV strain IbAr 2000 via the IP, intramuscular 
(IM), intranasal (IN) and SC routes, and showed that all animals 
died with an average time to death of 4+0, 5.20.6, 7x0 and 
4.640.2 dpi, respectively  (Zivcec et al., 2013) 
Thrombocytopenia,  coagulopathy, strong  pro-inflammatory 
responses were observed in these animals. Live CCHF of up to 
10^ TCIDso/mg of tissue could be detected in the blood, lymph 
node and various major organs (Zivcec et al, 2013). In 
contrast, infection of wild-type C57BL/6 mice with CCHFV did 
not result in any pathology (Zivcec et al., 2013). 


Severe fever with thrombocytopenia syndrome virus (SFTSV) 
Severe fever with thrombocytopenia is a newly recognized 
disease in rural areas of northeastern and central China, with 
several cases in Japan and South Korea (Promedmail, 2013). 
Caused by SFTSV, the transmission route of the virus is still 
unknown, but most likely involves arthropod vectors or animal 
hosts since the virus has been detected in ticks collected from 
domestic animals (Tian et al., 2017), and the animals (i.e., 
goats, cattle and dogs) also have high levels of SFTSV-specific 
antibodies (Jiao et al., 2012). Patients infected with SFTSV 
present with fever, vomiting, diarrhea, thrombocytopenia, 
leucopenia, and increased liver enzyme levels, in which severe 
cases of SFTSV eventually result in multiple organ failure 
resulting in death (Yu et al., 2011). The fatality rate amongst 
hospitalized patients can be up to 3096, and hundreds of cases 
are reported annually in China (Liu et al., 2015). 

Infection of wild-type mice (BALB/c, C57BL/6) results in 
limited weight loss but the animals do not succumb to disease 
(Chen et al., 2012; Jin et al., 2012). In one study, Liu et al. 
(2014) infected 6-10 week old /fnar^ mice (129/Sv 
background) SC with 109 ffu of SFTSV strain YL-1. The mice 
were highly susceptible to challenge, with all mice appearing ill 
by 3 dpi, resulting in death between 3-4 dpi. Blood and major 
organs (brain, heart, kidney, intestine, liver, lung and spleen) 
were collected from infected /fnar^ mice daily, and results 
showed high levels of virus replication with systemic spread to 
all organs. In particular, the spleen and intestine had the highest 
peak virus titers at death (Liu et al., 2014). In another study, 
Matsuno et al. infected 6-12 week old Ifnar^ mice (C57BL/6 
background) with either a high dose (10° TCIDso per animal) or 
a low dose (10? TCIDso per animal) of SFTSV strain SD4 via the 


intradermal (ID), IP, IM or SC routes. The results showed that 
the Ifnar^ mice were susceptible to infection via all routes, with 
animals succumbing to death at 4 and 6 dpi in the high and low 
dose groups, respectively (Matsuno et al., 2017). 


HENIPAVIRIDAE 


Hendra virus (HeV) and Nipah virus (NiV) 
HeV was discovered in 1994 as the etiologic agent that caused 
an acute respiratory disease in horses in Australia with sporadic 
but lethal transmission to humans, with one fatal case 
developing pneumonitis, respiratory and renal failure, arterial 
thrombosis, and eventually cardiac arrest seven days after 
admission (Selvey et al., 1995). HeV currently still poses a 
threat to Australian livestock, and the CFR is estimated to be 
60% for humans and 75% for horses (Field et al., 2011). NiV 
was discovered in 1999 in Malaysia with spread to neighbouring 
Singapore, resulting in 100 deaths from 257 human cases 
(CDC, 1999a) Patients typically present with respiratory 
problems and fever, as well as encephalitis with symptoms of 
headache, drowsiness, disorientation and confusion, rapidly 
progressing to coma. Since then, outbreaks of NiV have caused 
severe encephalitis in Bangladesh and India, with a CFR of 
~75% (Lo & Rota, 2008). Pigs are susceptible to infection and 
act as amplifying hosts to humans (CDC, 1999b). Fruit bats are 
the natural reservoir for both viruses (Halpin et al., 2011). 
Wild-type mice are only susceptible to HeV or NiV infection if 
the virus is administered via the intracranial (IC) route, but not 
through any other types of inoculations (Dhondt et al., 2013). In 
their study, Dhondt et al. infected 3-18 week old /fnar^ mice 
(C57BL/6 background) IP with 109 pfu of HeV. It was observed 
that while the infection was fully lethal in 3-week old mice, the 
susceptibility decreased with increasing age and the same dose 
of HeV in 18-week old mice only resulted in 5096 mortality. The 
moribund mice died between 7—13 dpi. For NiV strain UMMC1, 
4—12 week old /fnar^ mice (C57BL/6 background) were 
infected IP with increasing dosages from 100—109 pfu. The mice 
were found to be uniformly susceptible with deaths between 6— 
9 dpi in the 109 pfu group, and the LDso was calculated to be 
8x10? pfu in Ifnar^ mice. Infected mice with both viruses first 
showed behavioural changes including agitation, edginess and 
no grooming. Neurological symptoms were observed with 
advanced disease including tilted head and paralysis. A weight 
loss of approximately 1596—2596 was observed 1-2 days before 
death and found to be a good predictor of mortality (Dhondt et 
al., 2013). 


TOGAVIRIDAE 


Venezuelan equine encephalitis virus (VEEV) 

Venezuelan equine encephalitis was first identified in Venezuela 
in 1938, and outbreaks of the causative agent, VEEV, have 
occurred mostly in Central and South America, but the United 
States have also reported cases (Weaver et al., 2004). A 
mosquito-borne virus (Beaman & Turell, 1991), VEEV can infect 
and amplify in equine species, resulting in encephalitis as well 
as progressive disorders in the central nervous system. 
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Transmission of the virus to humans via the mosquito vector 
can result in the patient presenting with malaise, fever, 
headache and encephalitis (Weaver et al., 2004). The CFR is 
estimated to be 0.796—196 (Weaver et al., 1996), but permanent 
neurological damage have been noted with survivors of VEEV- 
induced encephalitis (León et al., 1975). 

Infection of wild-type mice with a virulent VEEV strain (V3000) 
results in death at approximately 8.3+0.5 dpi, but infection with 
an attenuated VEEV strain (V3032) does not result in mortality 
(Schoneboom et al., 2000). In their study, Schoneboom et al. 
infected 8-12 week old /fnar^ mice (129Sv/Ev background) 
with 1x10? pfu of VEEV either a virulent (V3000) or attenuated 
strain (V3032) SC into the left rear footpad. Within 18—20 hours 
of infection with either virus, the mice displayed hunching, 
ruffled fur and appeared lethargic. Advanced disease includes 
convulsions and prostration resulting in death. The mean time 
to death was 1 dpi for both VEEV strains, and high levels of live 
VEEV could be detected in the sera and brains of moribund 
Ifnar^ mice (Schoneboom et al., 2000). 


Chikungunya virus (CHIKV) 

First isolated in Tanzania in 1952 and carried by Aedes 
mosquitoes, CHIKV infections in humans result in Chikungunya 
fever, which is a severe illness in humans characterized by 
fever, headache, myalgia, rash, and acute as well as persistent 
arthralgia (Burt et al., 2012). Despite considerable morbidity, the 
CFR is estimated to be 0.196, with those older than 65 and/or 
underlying medical problems to be most at risk of death 
(Caglioti et al., 2013). CHIKV outbreaks occurred mostly in 
central/southern Africa and southeast Asia during the 1960s— 
2000s (Powers & Logue, 2007), but in recent years large-scale 
outbreaks have been reported on the island of Reunion (Roth et 
al., 2014), India (Pialoux et al., 2007) and also the Americas 
(Staples & Fischer, 2014), in which millions of infections were 
reported. 

Infection of wild-type adult mice with 10° pfu of CHIKV does 
not result in mortality (Couderc et al., 2008). In their study, 
Couderc et al. infected adult Ifnar~ mice (129s/v) via ID with 
20 pfu of CHIKV (Couderc et al., 2008) and found that the mice 
died within an average of 340.2 dpi, with an LDso of 3 pfu. 
Infectious virus could be detected in the liver within 16 hours 
after infection and abundantly detected in the muscles, joints, 
skin, brain, liver, spleen and sera by 3 dpi. Another study by Pal 
et al. infected 6-8 week old /fnar~ mice (C57BL/6 background) 
SC in the footpad with 20 ffu of CHIKV, and found that all mice 
died within 4 dpi, but that these animals did not develop the 
arthritis observed in humans (Pal et al., 2013). 


SUMMARY 


Immunocompromised /fnar^ mice have been shown to be a 
very good alternative small animal model for highly virulent 
pathogens that do not cause disease in immunocompetent 
mice. In this review, we described the different parameters and 
results from experimental infection of Ifnar^ mice with various 
pathogens (Table 1). It is obvious that the advent of /fnar1~ 
mice undoubtedly constituted a major step forwards in allowing 
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researchers to easily and rapidly study the pathogenesis of 
clinical isolates during a potential outbreak situation, as these 
animals are more susceptible to viral infections (Table 2). 

However, Ifnar^ mice as a model do have some weaknesses. 
Since these mice have defective innate immune responses 
which lead to impaired adaptive immunity, they are not good 
animal models for studying antiviral compounds, particularly 
vaccines (Züst et al., 2014). Additionally, it appears that the age 
of the /fnar* mice plays a role in host susceptibility to some 
viruses, as ZIKV loses the ability to cause disease in mice older 
than six weeks, whereas only partial lethality could be achieved 
with HeV infection in mice older than six weeks (Table 2). With 
respect to viruses from the Coronaviridae family, the /fnar~ 
mutation was found instead to prevent the lethal pneumonia 
observed in SARS-CoV mice (Channappanavar et al., 2016), 
whereas the absence of the human CD26 (a.k.a. DPP4) 
receptor for Middle East respiratory syndrome coronavirus 
(MERS-CoV) means that /fnar^- mice must be first transduced 
with a human adenovirus serotype 5-vector expressing human 
CD26 in order to become transiently susceptible to MERS-CoV 
infection via the IN route (Zhao et al., 2014). The transduced 
Ifnar^- mice were shown to experience ~20% body weight loss 
and delayed virus clearance by approximately 3 days compared 
to transduced wild-type mice, but do not die from the infection 
(Zhao et al., 2014). 

Aside from the transduction strategy, a popular method is to 
generate host-adapted viruses by sequential passaging in the 
livers and spleens of rodents (i.e., mice and guinea pigs) in 
vivo, in order to generate increasingly pathogenic virus variants 
that cause lethal disease to the immunocompetent host. This 
method has been widely used in the Filoviridae field to generate 
adapted viruses for EBOV (Bray et al., 1998; Volchkov et al., 
2000), MARV (Qiu et al., 2014), RAVV (Warfield et al., 2009) 
and SUDV (Wong et al., 2016) in wild-type mice or guinea pigs. 
In many cases, these viruses harbour very few mutations 
compared with the original clinical isolates, and the ability of the 
adapted virus to evade the host Type | IFN response (via 
mutations in the viral antigen responsible for this function) is 
positively correlated with its virulence in the host (Ebihara et al., 
2006). 

While animal models for studying virus pathogenesis leading 
to severe disease or lethality should always be the primary 
priority, an important future aim would be to also establish small 
animal models for studying pathogen transmission (initial work 
includes a guinea pig-based model of EBOV transmission 
(Wong et al., 2015), as well as developing small animal models 
to study various important phenomena of disease, such as the 
persistence of ZIKV in the testes of immunocompromised and 
immunocompetent mice (Govero et al., 2016; Ma et al., 2016), 
ZIKV infections leading to birth defects in wild-type mice 
(Cugola et al., 2016), ZIKV infections leading to microcephaly in 
neonatal mice (Li et al., 2016), or the adulthood sequelae of 
mice who survived congenital ZIKV infections (Cui et al., 2017). 
These studies in small animals will set the stage and provide 
important directives in subsequent investigations of similar 
disease phenomenon/sequelae in larger animal models, and 
ultimately, humans. 


Table 1 Summary of experimental parameters and results of Ifnar^ mice challenged with various outbreak viruses 





Pathogen Strain Age (weeks) and Challenge Challenge Death Mean time to death, References 
background dose route rate or range 
WNV 3000.0259 8-10, 129Sv/Ev 10? pfu SC 10096 Not provided Samuel & Diamond, 2005 
101 pfu Not provided 
10? pfu 3.8+0.5 days 
8-10, C57BL/6 10? pfu SC 10096  3.4+0.5 days 
DENV-2  PL046 5-6, 129Sv/Ev 10? pfu IV 096 N/A Shresta et al., 2004 
DENV-1  Mochizuki 4.4x10^ pfu N/A 
YFV Asibi 3-4, 129 10^ pfu SC 100% 7-9 days Meier et al., 2009 
Angola73 7-8 days 
JEV JEV JaOArS982 5-6, 129 10° pfu SC 9096 120 hours Aoki et al., 2014 
10? pfu 10096 96 hours 
10^ pfu 80 hours 
108 pfu 64 hours 
ZIKV MP1751 5-6, 129Sv/Ev 108 pfu SC 10096 6 days Dowall et al., 2016 
H/PF/2013 5-6, C57BL/6 10? pfu SC 8-10 days Lazear et al., 2016 
MR766 9—13 days 
Dakar 41671 10? ffu 6 days 
Dakar 41667 6 days 
Dakar 41519 6 days 
ZIKV FSS13025 3, C57BL/6 1x10? pfu SC 100% 6-7 days Rossi et al., 2016 
5, C57BL/6 50% 8-9 days 
11, C57BL/6 0% N/A 
EBOV Mayinga 8-16, 129 10° pfu SC 100% 7.3 days Bray, 2001 
IP 5.4 days 
Kikwit 0% N/A 
SUDV Boneface 100% 6.3 days 
RESTV 0% N/A 
TAFV 0% N/A 
MARV Musoke 67% 8.5 days 
RAVV 100% 6.0 days 
MARV Popp 6-9, 129 1028 TCIDso Aerosol 100% 13.0 days Lever et al., 2012 
1038 TCIDso 12.0 days 
1048 TCIDso 10.2 days 
1098 TCIDso 11.0 days 
EBOV E719 10? TCIDso 8.0 days 
10! TCIDso 8.0 days 
10? TCIDso 8.0 days 
SUDV Boneface 104° TCIDso 096 N/A 
LASV Josiah 8-12, 129/Sv 10? ffu IV 096 N/A Rieger et al., 2013 
AV 
BA366 
Nig04-10 
CCHFV . IbAr2000 7—10, 129 Sv/Ev 10! ffu IP 10096 4 days Bereczky et al., 2010 
10? ffu 3 days 
10? ffu 2 days 
108 ffu 2 days 
6-12, C57BL/6 104 TCIDso IP 100% 4+0 days Zivcec et al., 2013 
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Continued 





Pathogen Strain Age (weeks) and Challenge Challenge Death Mean time to death, ON. 
background dose route rate or range 
CCHFV IbAr2000 6-12, C57BL/6 10^ TCIDso IM 100%  5.2x0.6 days Zivcec et al., 2013 
IN 7+0 days 
SC 4.6+0.2 days 
SFTSV YL-1 6-10, 129/Sv 108 ffu sc 100% 3-4 days Liu et al., 2014 
SD4 6-12, C57BL/6 10? TCIDso IP 100% 5 days Matsuno et al., 2017 
IM 5-6 days 
SC 5 days 
ID 6 days 
105 TCIDso — IP 3-4 days 
IM 4 days 
SC 4 days 
ID 6 days 
HeV SD4 3, C57BL/6 108 pfu IP 100% 11 days Dhondt et al., 2013 
6, C57BL/6 83% 11-13 days 
18, C57BL/6 50% 7 days 
NiV UMMC1 4—12, C57BL/6 10? pfu 096 N/A 
10? pfu 1796 10 days 
104 pfu 67% 10 days 
10° pfu 83% 8-10 days 
108 pfu 100% 6-9 days 
VEEV V3000 8-12, 129Sv/Ev 1x10? pfu SC 100% 1 day Schoneboom et al., 2000 
V3032 1 day 
CHIKV 21 Adult (age not given), — 20 pfu ID 100% 320.2 days Couderc et al., 2008 
129s/v 
LR 6-8, C57BL/6 20 ffu SC 100% 3-4 days Pal et al., 2013 


IM: Intramuscular; IN: Intranasal; SC: Subcutaneous; IP: Intraperitoneal; ID: Intradermal; IV: Intravenous. 


Table 2 Advantages and disadvantages of using /fnar^- mice for studying human infectious diseases, compared to other strategies 


Strategy in small animal models 


Advantages 


Disadvantages 





Knockout mice 


Virus adaptation to host via 
sequential passaging 


Transduction with adenoviral 
vectors encoding the entry 
receptor to confer sensitivity 


Susceptible to a wide range of clinical 
isolates of viruses 

Can study pathogenesis of a new pathogen 
rapidly 


Can cause uniform lethality 
Good for screening drugs and vaccines 


Wild-type mice are widely available 


Useful when no other known small animal 
models exist (i.e., MERS) 
Can test with clinical isolate of virus 


MERS: Middle East respiratory syndrome. 
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Cannot study immune responses properly due to abnormal 
innate immunity 
Some viruses may not cause disease in knockout mice 


Cannot test drugs and vaccines effectively 


Can be age-sensitive: older mice may lose their susceptibility 
to the pathogen 
Not always successful in creating a lethal variant 


Can be time consuming to create a lethal variant 


Not clinical isolate of virus and thus may harbour important 
differences in pathogenesis 
Need to know the identity of the receptor 


Time consuming to create the recombinant adenovirus 
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ABSTRACT 


The family Filoviridae, which includes the genera 
Marburgvirus and Ebolavirus, contains some of the 
most pathogenic viruses in humans and non-human 
primates (NHPs), causing severe hemorrhagic 
fevers with high fatality rates. Small animal models 
against filoviruses using mice, guinea pigs, 
hamsters, and ferrets have been developed with the 
goal of screening candidate vaccines and antivirals, 
before testing in the gold standard NHP models. In 
this review, we summarize the different animal 
models used to understand filovirus pathogenesis, 
and discuss the advantages and disadvantages of 
each model with respect to filovirus disease 
research. 


Keywords: Filovirus; Ebola virus; Marburg virus; 
Marburg virus disease; Ebola virus disease; Mouse- 
adapted ebolavirus; Guinea-pig-adapted ebolavirus; 
Sudan virus; Plaque-forming units 


INTRODUCTION 


The family Filoviridae consists of three genera, Ebolavirus, 
Marburgvirus, and Cuevavirus (Kuhn et al. 2011, 2013). 
Filoviruses are classified as biosafety level 4 (BSL-4) 
pathogens and infections result in severe hemorrhagic fevers in 
humans and non-human primates (NHPs), with fatality rates as 
high as 9096 (Sanchez et al., 2007). Marburgvirus has two 
member viruses: Marburg virus (MARV) and Ravn virus (RAVV). 
Ebolavirus has five member viruses: Ebola virus (EBOV), 
Sudan virus (SUDV), Tai Forest virus (TAFV), Bundibugyo virus 
(BDBV), and Reston virus (RESTV) (Kuhn et al., 2013). Live 
viruses belonging to Cuevavirus have not yet been isolated. 
Currently, no approved vaccines or therapeutics are available 
against filovirus infections; however, clinical trials of candidate 
vaccines and therapeutics conducted during the 2014—2016 
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EBOV outbreak in West Africa (Martin et al., 2016; Wong & 
Kobinger, 2015) provide hope that a licensed medical 
countermeasure is on the horizon. 

The transmission of filoviruses to humans is poorly 
understood, but outbreaks are likely started by direct contact 
with dead animal carcasses (typically NHPs) or in the case of 
MARV, the bat reservoir, resulting in subsequent transmission 
within the susceptible human population (Leroy, 2009). The 
blood and/or other body fluids of dead animals and humans are 
especially infectious, with viremia in human patients reaching 
10? copies/mL of blood (Faye et al., 2015). Contact with blood 
and bodily fluids is thought to be the main mode of 
transmission. After an incubation period of 2-21 days, general 
symptoms such as fever, chills, fatigue, headache, and 
myalgia appear. As disease progresses, systemic (prostration, 
lethargy), gastrointestinal (anorexia, vomiting, abdominal pain, 
diarrhea), respiratory (chest pain, breath shortness, cough, 
nasal discharge) vascular (conjunctival injection, postural 
hypotension, edema), and neurological (headache, confusion, 
coma) symptoms may appear, sometimes accompanied with 
hemorrhage from venipuncture sites. In fatal cases, patients 
succumb after multiple organ failure between 6-16 days after 
the onset of symptoms (Nakayama & Saijo, 2013). 

Animal models against filoviruses have been developed in 
mice, guinea pigs, hamsters, ferrets and NHPs (Bente et al., 
2009; Bradfute et al., 2012; Connolly et al., 1999; Marzi et al., 
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2016; Kozak et al., 2016; Kroeker et al., 2017b; Yamaoka et al., 
2017). The development of animal models that recapitulates 
hallmarks of human filoviral disease is crucial in understanding 
the pathogenesis of these viruses. In this review, we summarize 


Table 1 Animal models for studying filovirus infections 


the animal models that have been developed for filoviruses 
(Table 1), and discuss the advantages and disadvantages of 
each animal model (Table 2). 





Filovirus species Immuno-competent mouse Immuno-compromised mouse Guinea pig Syrian hamster Ferret NHP 
EBOV + + + + + + 
MARV + + + + - + 
RAVV + + + - - + 
SUDV - + + - + + 
TAFV - - - - - - 
BDBV - - - - + + 


NHP: Non-human primate; EBOV: Ebola virus; MARV: Marburg virus; RAVV: Ravn virus; SUDV: Sudan virus; TAFV: Tai Forest virus; BDBV: 


Bundibugyo virus; “+”: Available; "-": Not available. 


Table 2 Advantages and disadvantages of animal models for filovirus infections 


Disadvantages 





Animal models Advantages 
Mice Low cost, easy to use 

Transgenic and knockout models are available 
Guinea pigs Low cost, larger animals to study disease 


progression and easy to use 


Syrian Hamsters — An alternative to guinea pigs 
Elegant model to compare differences in immune 


responses to EBOV infection 


Ferrets Small body size, low cost 
Closely resembles to human filoviral disease 
Can use clinical isolates to study pathogenesis 
NHPs Gold standard model to evaluate filovirus 


infections and closely recapitulate human disease 


NHPs: Non-human primates. 


MOUSE MODELS 


Mice are easy to handle in the laboratory, are available 
commercially in large numbers, and have a low unit cost. Due 
to the wide availability of biochemical reagents and immunological 
tools, mice are preferred as the small animal model for filovirus 
research (Bradfute et al., 2007). 


Immunocompetent Mouse Model 

Ebola virus (EBOV) 

As immunocompetent mice are resistant to infections with wild 
type EBOV (WT-EBOV) (Banadyga et al., 2016; Bray, 2001; 
Warfield et al., 2009), Bray et al. (1998) adapted the 1976 
Mayinga isolate of EBOV to cause lethal disease in adult 
immunocompetent mice through the sequential passaging of 
infected liver homogenates in newborn mice. The resultant 
plaque-purified mouse-adapted EBOV (MA-EBOV) was lethal 
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Only i.p. infection is 100% lethal 
Mouse-adapted variants needed 


Transgenic and knockout models are not available 

Lack of immunological tools and reagents to evaluate cell-mediated 
responses to vaccines and therapeutics 

Guinea pig-adapted variants needed 


Mouse-adapted variants needed 
Lack of commercially available reagents and host genomic information 
Transgenic and knockout models are available 


Limited availability of ferret-specific reagents 
Ferret immune responses are poorly understood 
Transgenic and knockout models are available 


Animals are expensive, ethical considerations and extensively 
husbandry requirements needed 
Transgenic and knockout models are not available 


for adult BALB/c, CD-1, and C57BL/6 mice at 30 or 3 000 
median lethal dose (LDso) of 10° or 10? plaque-forming units 
(pfu) when injected intraperitoneally (i.p.). Remarkably, infection 
with MA-EBOV did not produce any disease symptoms via the 
subcutaneous (s.c.) or intramuscular (i.m.) routes at doses of 
10° pfu (Bray et al., 1998). Compared with WT-EBOV, MA- 
EBOV exhibited eight amino acid changes in the coding and 
non-coding regions of the virus genome (Ebihara et al., 2006), 
with amino acid changes observed in viral protein (VP) 35, 
VP24, nucleoprotein (NP), and RNA-dependent RNA polymerase 
(L) genes. Furthermore, the determinants of virulence in mice 
were shown to be mutations in NP and VP24, in which 
recombinant viruses containing only the NP and VP24 
mutations were resistant to type | interferon (IFN) in vitro and 
conferred lethality to mice (Ebihara et al., 2006). The disease 
observed in adult immunocompetent mice by i.p. inoculation of 
MA-EBOV closely reflected EBOV infection in humans, and the 
virus replicated well in blood, reaching peak titers of 10° pfu/mL. 


Biochemical and  histopathological studies have shown 
diminished functions in the liver and kidney (Gibb et al., 2001b). 
In addition, lymphocyte apoptosis and decrease in platelet 
counts were observed in mice infected with MA-EBOV. In livers 
of MA-EBOV-infected mice, viral replication was seen in 
hepatocytes, Kupffer cells, and sinusoidal endothelial lining 
cells. In the spleen, viral antigen was noted at day 2 after 
infection. Coagulopathy, such as disseminated intravascular 
coagulation with prolongation of prothrombin time (PT) and 
activated partial thromboplastin time (APTT) have also been 
detected at advanced stages of the disease (Bray et al., 2001; 
Warfield et al., 2009). However, fibrin deposition and 
breakdown at sites of viral replication in the spleen and other 
tissues have not been observed in tissue sections of MA- 
EBOV-infected mice (Bray et al., 1998). The MA-EBOV initially 
targets macrophages and other mononuclear phagocytes to 
invade the regional lymph nodes, similar to observations in 
humans, NHPs, and guinea pigs (Connolly et al., 1999; Davis et 
al., 1997; Gibb et al., 2001a, b; Zaki & Goldsmith, 1999). Thus, 
immunocompetent mouse models are considered ideal for 
studying fundamental aspects of viral replication, pathogenesis, 
and counter immune responses. They also provide an excellent 
platform for evaluating the therapeutic efficacy of candidate 
vaccines, antibodies, and small molecules such as anti-viral 
drugs (Kroeker et al., 2017a; Mendoza et al., 2016; Qiu et al., 
2016; Zeitlin et al., 2016). 


Marburg virus (MARV) 

Marburg virus (MARV) and RAVV variants lethal to adult 
immunocompetent mice have been developed through 
sequential passage in the livers and spleens of severe, 
combined immunodeficient (SCID) mice (Warfield et al., 2007). 
Qiu et al. (2014) established a lethal mouse model against the 
Angolan strain MARV (MARV/Ang), the most virulent strain of 
MARY, with mice infected with 2 000 LDso of MARV/Ang-MA i.p. 
showing lymphopenia, thrombocytopenia, and uncontrolled 
viremia. Warfield et al. (2009) developed a lethal mouse- 
adapted variant of RAVV (MA-RAVV), which caused death 5— 
10 days post infection (dpi), with disease manifestations 
including high viremia, fibrin degradation products, platelet loss, 
profound loss of tissue lymphocytes, and liver damage 


Immunocompromised Mouse Model 
Ebola virus (EBOV) 
Another method of studying WT filovirus infections in mice 
involves the use of mice with defects in the immune system. 
Bray et al. (2001) reported that i.p. inoculation of WT-EBOV 
was lethal to IFN-a/B receptor knockout (IFNAR"^) mice, with 
immunohistochemical analysis demonstrating increased 
replicating virus and viral antigen in livers and spleens. 
Additionally, mice knocked out for transcription factor STAT1 
(which plays a key role in type | IFN signaling) were also highly 
susceptible to i.p. inoculation with 10? pfu of WT-EBOV (Bray et 
al., 2001). 

The type | IFN response plays a key role to resistance in 
mice. In addition, mouse models deficient in type | IFN signaling 
are more susceptible to WT filoviruses, serving an important 


model for studying pathogenesis. Immunodeficient mouse 
models of EBOV infection demonstrate high viremia and viral 
load in the spleen, liver, and multiple organ tissues. 
Lymphopenia, thrombocytopenia, kidney dysfunction, and liver 
damage have been observed in disease progression (Bradfute 
et al., 2007; Bradfute et al., 2010; Bray et al., 2001; Shurtleff & 
Bavari, 2015). As immunocompromised mice have a higher unit 
cost compared with immunocompetent mice, require handling 
in sterile conditions, and cannot mount a normal robust immune 
response (thus cannot determine immune correlates of 
protection or pathogenesis), they have not been used as 
extensively for research. However, the advantage to this animal 
model is that clinical WT filovirus strains can be used as a 
challenge virus to test the efficacy of candidate 
countermeasures rapidly without the need for host adaptation, 
which can be time consuming. 


Marburg virus (MARV) 

Warfield et al. (2007) examined SCID mice (which lack B and T 
cells) with intact IFN responses, who became sick when 
inoculated i.p. with ~10° pfu of ‘SCID mouse-adapted’ MARV- 
Musoke and died within 3-4 weeks after infection. However, 
after subsequent animal-to-animal passages, both RAVV and 
Musoke strains of MARV (MARV-Musoke and MARV-Ci67) 
caused rapidly lethal disease (Warfield et al., 2007). Initial signs 
of MARV disease in these animal models include fever, 
anorexia, rash, huddling, weight loss, dehydration, and 
diarrhea. Severe complications such as prostration, failure to 
respond to stimulation, hind limb paralysis, and bleeding from 
body orifices develop at 6—10 days. Early hematological and 
immunological changes include lymphopenia, neutrophilia, and 
profound thrombocytopenia, whereas alterations in serum 
chemistry levels, such as increases in liver enzymes, are also 
prominent after infection. At death, MARV was found at high 
titers in the blood, liver, spleen, kidneys, and other major 
organs, indicating systemic spread of the virus (Warfield et al., 
2007). 


Sudan virus (SUDV) 

Mice deficient for type 1 IFN infections (a/B/y IFN receptor 
knockout, IFN-a/BR^) have been found to be susceptible to 
wild-type SUDV (Brannan et al., 2015). Specifically, IFN-a/BR 
mice challenged i.p. with a dose of 10? pfu of SUDV Boneface 
variant (SUDV/Bon) subsequently exhibited ~25% weight loss, 
which peaked at 7-9 dpi. Maximum lethality was observed in 
the 10? pfu challenge group, with a 6396 survival, and mean time 
to death was approximately 11 days. Histological analyses of liver 
and spleen samples of the SUDV-infected mice on day 3 showed 
the development of hepatocellular degeneration and necrosis. 


GUINEA PIG MODEL 


Ebola virus (EBOV) 

Guinea pigs occupy a prominent role in filovirus research and 
are widely used for studying filoviral hemorrhagic fevers (Bowen 
et al., 1977; Cross et al., 2015a; Wong et al., 2015a, b, c). 
However, guinea pigs infected with WT filoviruses result in only 
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transient illness (Ryabchikova et al., 1996; Simpson et al., 
1968). Therefore, EBOV was serially passaged in the livers and 
spleens of guinea pigs until the adapted virus (GPA-EBOV) 
caused uniform lethality, with animals succumbing to death at 
7—9 dpi (Simpson et al., 1968). At a dose of 103.8 pfu, GPA- 
EBOV-infected animals showed fibrin deposition and 
thrombocytopenia during the late stages of infection (Connolly 
et al., 1999), as well as high viral titers in the spleen, liver, 
adrenal gland, and lungs. Viremia was observed in guinea pigs 
2 days after inoculation, which peaked on day 7 (>10* pfu/mL) 
(Connolly et al., 1999; Subbotina et al., 2010). 

Guinea pigs infected with GPA-EBOV showed drastic 
histopathological changes in Kupffer cells (especially virus 
replication), death of hepatocytes, and destruction of the 
sinusoid wall (Connolly et al., 1999), with the spleen and lymph 
nodes demonstrating lymphoid necrosis along with neutrophilia 
and lymphopenia (Connolly et al., 1999; Subbotina et al., 2010). 
Remarkably, lymphocyte bystander apoptosis was not 
prominent in these animals (Bradfute et al., 2007; Bray et al., 
1998; Connolly et al., 1999). While the guinea pig model 
recapitulates some aspects of human filoviral disease, its higher 
unit cost, need for more specialized housing/husbandry 
requirements, and the lack of reagents to characterize aspects 
of innate/adaptive immune responses (such as T-cell immunity 
or cytokine levels) mean that it is more suited as a secondary 
animal model for confirming experimental results and trends 
from mice studies. 

Adaptation of MARV (isolates Ci67 and Musoke) in guinea 
pigs from previous research resulted in their deaths on days 7— 
9, with body temperatures recorded at 41.1 °C (Simpson et al., 
1968) and clinical symptoms such as bloated face and loss of 
appetite and weight evident. Cross et al. (2015b) evaluated and 
compared pathogenicity in outbred guinea pigs with GPA- 
RAVVutme and GPA-MARV infections to better understand the 
disease pathogenesis of MARV infection in the guinea pig 
model. Results demonstrated prolonged clotting times (PT and 
APTT) associated with an increase in plasminogen activator 
inhibitor 1 and von Willebrand factor activity in GPA-RAVV and 
GPA-MARV-Angola-infected animals. In addition, circulating 
protein C activity and tissue factor concentrations also 
decreased during the time of death. 


Marburg virus (MARV) 

Guinea pigs infected with approximately 5x10? pfu of GPA- 
MARV-Angola showed remarkable differences in tissue factor 
kinetics (Cross et al., 2015b). In these animals, prolonged 
clotting times (PT and APTT) were associated with an increase 
in plasminogen activator inhibitor 1 and von Willebrand factor 
activity. Circulating protein C activity and tissue factor 
concentrations were also decreased at death. Thus, these 
proteins likely played a significant role in coagulopathy. 
Furthermore, disturbance of the fibrinolysis pathway was 
implicated by the depression in the metabolism of bradykinin, 
prekallikrein, and thrombin-activated fibrinolysis inhibitor (TAFI). 
Interestingly, GPA-RAVV-infected animals displayed increased 
TAFI concentrations at 3 dpi, whereas GPA-MARV-Angola- 
infected animals showed a decrease in TAFI concentrations at 
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late stage infection. Inflammatory mediators such as TNF-a, IL- 
6, nitric oxide, and HMGB-1 were also elevated during the late 
stages of infection. 


Sudan virus (SUDV) 

Wong et al. (2015a) developed and characterized a guinea 
pig-adapted SUDV variant based on the Boneface variant 
through serial passaging of WT-SUDV in the livers and 
spleens, and showed that the adapted virus caused uniform 
lethality in guinea pigs. The guinea pig-adapted SUDV 
variant (SUDV-GA) caused LDso to animals at a dose of 
5.3x10% 50% tissue culture infective doses (TCIDs50). 
Animals infected with SUDV-GA showed high viremia and 
died between 9-14 dpi. Several hallmarks of SUDV infection 
were observed, such as lymphadenopathy, elevated liver 
enzyme activities, and coagulation abnormalities. 
Histopathology and immunohistochemistry findings indicated 
that the SUDV-GA antigen was present in the livers and 
spleens of infected animals. 


SYRIAN HAMSTER MODEL 


Ebola virus (EBOV) 

A Syrian golden hamster model for EBOV was developed by 
inoculating 6-week-old hamsters i.p. with 10? focus forming 
units (ffu) of MA-EBOV (Ebihara et al., 2013; Zivcec et al., 
2011). Hamsters i.p. inoculated with MA-EBOV (Bray et al., 
1998) showed signs of disease, such as ruffled fur and 
decreased activity beginning at 3 dpi, and all animals died by 
days 4—5, with correspondingly high virus titers in the spleen, 
liver, kidneys, heart, lungs, and brain, indicating systemic 
spread of the virus (Ebihara et al., 2013). 

Viremia was detected beginning on day 2 and reached -10? 
pfu/mL by day 4, with infected hamsters showing rapid and 
increased distribution of the viral antigen in the spleen and liver 
(Ebihara et al, 2013) In the spleen, macrophages and 
marginal reticular-like cells in the red pulp and marginal zone 
showed higher virus replication. In the liver, Kupffer cells were 
the first target cells for MA-EBOV, which then spread to 
hepatocytes starting at 2 dpi. Importantly, coagulopathy was 
observed in the MA-EBOV-infected Syrian hamsters, but was 
absent in other rodent models. In addition, thrombocytopenia 
and decreased protein C concentrations were also noted in the 
infected hamsters. Using in-house primer-probe sets, Ebihara 
et al. (2013) also showed that MA-EBOV infection induced the 
expression of cytokine/chemokine genes, including IL-1b, IL-6, 
TNF-a, IL-12p35, IP-10, and IL-10, in the spleen and liver, 
indicating an uncontrolled immune response. Due to the 
presence of rash and induction of cytokines/chemokines, the 
Syrian hamster more closely recapitulates EBOV disease 
compared with guinea pigs; however, the higher unit costs, 
need for specialized housing/husbandry equipment, and lack of 
commercial reagents for studying innate and cell-mediated 
immune responses means that this model is still not widely 
used in filovirus research. 


Marburg virus (MARV) 

Marzi et al. (2016) recently developed a novel Syrian golden 
hamster model for MARV infection using a hamster-adapted 
MARV variant Angola (HA-MARV). This model provides new 
insight on virus pathogenesis in hamsters, and moreover 
closely resembles MARV infection in humans and NHPs, 
including hemorrhagic manifestations and coagulation 
abnormalities. In this study, hamsters were infected i.p. with 10-— 
10? pfu of HA-MARV, leading to significant weight loss and 
disease surrender between 7—11 dpi. The HA-MARV infected 
hamsters displayed an increase in body temperature at 6 dpi 
and drop in temperature during the end stage of the disease. In 
contrast, the mock-infected animals (WT-MARV) did not show 
any changes in body temperature or weight, suggesting that 
infection with WT-MARV failed to induce the disease in 
hamsters. Maculopapular rashes along with visible petechial 
hemorrhages on the face and abdomen were noted in the HA- 
MARV-infected hamsters on 7 dpi. Hemorrhage was also 
identified in the footpads and joints, as well as in the kidneys. 
Coagulation parameters in HA-MARV-infected hamsters such 
as PT, APTT, and thrombin increased with delay in blood clot 
formation during the infection. However, the animals infected 
with WT-MARV did not exhibit these coagulation abnormalities. 
The HA-MARV replication was found to be robust, with viremia 
peaking at 108 TCIDso/mL in the blood, 107 TCIDso/mL in the 
liver, 105 TClDso/mL in the spleen, and 10° TCIDso/mL in the 
mesenteric lymph nodes by 7-8 dpi. Dysregulated immune 
responses were also observed in hamsters infected with HA- 
MARV; for example, pro-inflammatory chemokines MIP-1a and 
IP-10, as well as type | interferon responses, were elevated 
when compared to WT-MARV-infected animals. 


FERRET MODEL 


Ebola virus (EBOV) 

The susceptibility of domestic ferrets (Mustela putorius furo) to 
EBOV has been studied recently. Cross et al. (2016) 
demonstrated that ferrets infected intranasally (i.n.) with the 
EBOV Kikwit variant at 10° pfu, showed uncontrolled viral 
replication, abnormalities in white blood cell counts, and 
multiple-organ failure. Initial signs of the disease in the EBOV- 
infected ferrets were first observed at 3 dpi, followed by rapid- 
onset hypothermia and weight loss on day 4. Remarkably, 
ferrets died between 6-9 dpi. Necropsy inspection revealed 
pathological features of hemorrhagic fever, including petechial 
rashes on the skin surface, reticulated pallor of the liver, and 
mottled splenomegaly. All ferrets showed progressive neutrophilia 
and lymphocytopenia beginning on 4 dpi. Vascular leakage was 
also noted with hypoalbuminemia and hypoproteinemia, and 
increasing levels of circulating proinflammatory markers TNF-a 
and nitric oxide were also recorded on day 4. 

Kozak et al. (2016) also reported on a ferret model for EBOV. 
Animals were infected with 200 TCIDso via the i.m or i.n route, 
and showed peak viremia of 10" GEQ/mL or 10? TCIDso/mL by 
5 or 6 dpi, with evidence of sporadic viral shedding from the 
oral, nasal, and rectal cavities. Notably, weight loss was less 
pronounced in animals infected with EBOV compared with 


BDBV-infected animals from the same study. Biochemical 
analyses further demonstrated increases in concentrations/ 
activities of ALT, ALP, and BIL, suggesting liver damage, and 
increases in BUN and CRE, indicating elevated renal problems. 
Prolonged APTT and thrombin time (TT), along with an 
increase in fibrinogen levels, were also noted, suggesting that 
ferrets recapitulated disseminated intravascular coagulation. 
Interestingly, EBOV-infected ferrets displayed decreased levels 
of serum albumin after EBOV challenge, which might be the 
cause of disease-induced edema. 

The results from the above studies demonstrate that ferrets 
are highly susceptible to EBOV infection and the pathological 
parameters observed during disease pathogenesis closely 
resemble those found in humans. 


Sudan virus (SUDV) 

Kroeker et al. (2017b) demonstrated that ferrets are susceptible 
to wild-type SUDV infections. Ferrets inoculated with WT-SUDV 
via i.m and i.n experienced viremia, organ dysfunction, viral 
shedding, and eventual death. A decrease in white blood cell 
counts, increase in fibrinogen, and decrease in platelets and 
PT%, which are indicative symptoms of disseminated 
intravascular coagulopathy, were observed further during 
disease progression. Cross et al. (2016) also described a lethal 
model of SUDV with the domestic ferret. Ferrets were 
challenged i.n. with 10? pfu of SUDV, with initial signs of the 
disease observed at 3 dpi and rapid-onset hypothermia and weight 
loss detected by day 4 for SUDV-infected ferrets. Clinical signs 
included diarrhea, dehydration, nasal and ocular discharge, 
labored breathing, hunched posture, and altered gait, with 
death occurring at 6—9 dpi. Hematoxylin-eosin staining revealed 
that the most significant lesions in the infected ferrets were 
lymphohistiocytic and neutrophilic necrotizing hepatitis and 
necrotizing splenitis. 


Bundibugyo virus (BDBV) 

Kozak et al. (2016) challenged ferrets with 159 TCIDso of BDBV 
via the i.m. route, and observed viremia at 4 dpi and uniform 
lethality at 8 dpi. In all BDBV-infected animals, significant 
increases in coagulation parameters were observed, including 
in APTT, TT, and PT. Hematological analysis revealed a 
decrease in white blood cell counts, lymphocytes, and platelets 
in ferrets after infection. Petechial rashes were also observed 
during the disease. 


NON-HUMAN PRIMATE (NHP) MODELS 


Ebola virus (EBOV) 

Many reports have demonstrated African green monkeys 
(Chlorocebus aethiops), cynomolgus macaques (Macaca 
fascicularis), and rhesus macaques (Macaca mulatta) as 
excellent models for EBOV infection (Bowen et al., 1978; Davis 
et al., 1997; Ellis et al., 1978; Fisher-Hoch et al., 1992; Geisbert 
et al., 2003a, b, c, d; Ignatiev et al., 2000; Jaax et al., 1996; 
Jahrling et al., 1996; Johnson et al., 1995; Ryabchikova et al., 
1999; Wong et al., 2016). Following EBOV infection, NHPs 
became febrile (above 40 ?C) at 3 dpi, with pyrexia and 
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temperature decrease exhibited throughout the course of the 
disease, followed by death within 5-8 dpi (Baskerville et al., 
1978; Bowen et al., 1978; Ellis et al., 1978; Fisher-Hoch et al., 
1985; Luchko et al., 1995). Viremia has been observed at 3 dpi 
with titers reaching 1095-10" pfu/mL on days 4—5 (Bowen et al., 
1978; Fisher-Hoch et al., 1992; Jahrling et al., 1996). Compared 
with other NHP species, cynomolgus and rhesus macaques are 
the best available animal models due to their susceptibility to 
EBOV infections (Bente et al., 2009), thus providing a stringent 
test for any candidate medical countermeasure. In addition, 
African green monkeys and baboons infected with EBOV show 
impairment of the coagulation system. For example, fibrin 
thrombosis was localized to all visceral organs in African green 
monkeys, whereas hemorrhage was most prominent in the 
visceral organs such as the liver and spleen in baboons 
(Ignatiev et al., 2000; Ryabchikova et al., 1999). 

Currently, macaques well recapitulate pathological features of 
fatal filovirus disease observed in humans, including high 
viremia, coagulation abnormalities such as decreased platelet 
levels and increased blood clotting times, and aberrant 
proinflammatory cytokine responses, including the release of IL- 
6 and MCP-1 (Geisbert et al., 2003a, 2015). Studies in 
macaques have also demonstrated that doses of EBOV as low 
as 2-15 pfu administered by different challenge routes can 
produce lethal filovirus infection (Sullivan et al., 2000, 2003). 
Ebola virus disease (EVD) signs in macaques initially occur at 
3-5 days after exposure, with symptoms including fever and 
malaise, followed by anorexia, depression, lethargy, diarrhea, 
vomiting, and maculopapular rash. Hemorrhagic manifestations 
have also been identified, and include petechiae, ecchymosis, 
and bruising, hemorrhage at venipuncture sites, epistaxis, 
hematochezia, and hematuria. In addition, neutrophilia, 
lymphopenia, thrombocytopenia, and early monocytosis 
indicate abnormalities in complete blood count parameters. 

Although airborne transmission is not considered a significant 
route of human infection (Osterholm et al., 2015), aerosolized 
viruses have caused lethal disease in experimentally-infected 
NHPs. Previous study showed that control rhesus macaques, 
located 3 m from experimental rhesus macaques inoculated 
i.m. with EBOV, became infected (Jaax et al., 1995), with 
antigen staining patterns of pulmonary specimens implicating 
aerosol infection. Conversely, other research has shown no 
transmission between rhesus macaques infected i.m. with 
EBOV and uninfected rhesus macaques when housed adjacent 
to each other (Alimonti et al., 2014). However, virus 
transmission from pigs infected oro-nasally with EBOV to naive 
NHPs without direct contact was shown to be possible, with the 
infected NHPs exhibiting extensive lung damage (Weingartl et 
al., 2012). Mucosal exposure as a source of infection in rhesus 
macaques (through conjunctival and oral routes), required higher 
doses (5.2 logio of EBOV Mayinga isolate) compared to 
parenteral routes; however, doses of 10 pfu by oral or 
conjunctival routes did not result in any clinical disease (Mire et 
al., 2013). Thus, despite the above studies, an aerosol model of 
EBOV infection in NHPs consistently causing 100% lethality has 
not yet been established. 
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Marburg virus (MARV) 

NHPs are susceptible to MARV infections and have been 
extensively used to study its pathogenesis. Fernando et al. 
(2015) reported that cynomolgus macaques showed severe 
disease against Marburg Angola compared with other MARV 
strains, with high fever, anorexia, and lymphopenia, followed by 
death (mean time after infection of 6.7 days). Peak viremia 
reached 104-10’ TCIDso/mL by 4 dpi. Cynomolgus macaques 
showed febrile illness, anorexia, diarrhea, skin rash, and 
hemorrhagic manifestations at 2-6 dpi following high dose (10? 
pfu) infection by MARV Angola (Alves et al., 2010; Geisbert et 
al., 2007; Hensley et al., 2011; Simpson et al., 1968; Simpson, 
1969; Murphy et al., 1971), with a sudden decrease in body 
temperature followed by death at 6—13 dpi. Lymphocytosis was 
observed at the beginning of the illness (Geisbert et al., 2007; 
Gonchar et al., 1991; Johnson et al., 1996; Simpson et al., 1968; 
Simpson, 1969; Spiridonov et al., 1992), whereas thrombocytopenia 
and leukocytosis with increased neutrophilia were observed 5— 
6 dpi (Hensley et al., 2011). The infection was also observed to 
spread to the liver, adrenal glands, and finally to endothelial 
cells in a variety of organ tissues (Hensley et al., 2011). Viremia 
occurred on day 3 in cynomolgus macaques and African green 
monkeys, reaching peak titers of 107-108 pfu/mL at 8 dpi 
(Hensley et al., 2011). 


Sudan virus (SUDV) 

SUDV infection has been shown to cause death in the NHP 
model (rhesus macaques) 7-10 dpi (Thi et al., 2016). 
Specifically, adult rhesus macaques inoculated i.m. with 10° pfu 
of SUDV Gulu succumbed to infection after 7-10 days, with 
viremia peaking at 109? pfu/mL at 7 dpi, accompanied by liver 
and renal dysfunction as evidenced by increased ALT, GGT, 
BUN, and CRE levels. Control and treated animals that 
succumbed to SUDV infection showed lesions (observed by 
hematoxylin and eosin staining) in tissues. Significant lesions 
included splenic lymphoid depletion with tangible body 
macrophages and fragmented nuclear debris and enlargement 
of splenic red pulp with fibrin, multifocal necrotizing hepatitis 
with sinusoidal leukocytosis, and mild interstitial pneumonia. 
Furthermore, depression, anorexia, petechial rash, and 
hemorrhage were observed in the animals infected with SUDV. 


Bundibugyo virus (BDBV) 

Cynomolgus macaques have been used for studying the 
pathogenesis of BDBV infection (Mire et al., 2013), in which 
macaques were challenged i.m. with 10° pfu of BDBV. Animals 
were monitored for clinical signs of illness over the course of 28 
days. Clinical scores were recorded each day post-challenge 
for each animal using a scoring system based on dyspnea, 
depression, recumbency, and rash. No antibody titers (IgG) for 
the control group were observed against BDBV. The signs of 
disease in response to BDBV infection were more dramatic for 
the control animals. Signs observed between day 0-28 after 
BDBV challenge included fever, anorexia, depression, 
rectorrhagia, lymphopenia, and thrombocytopenia. Serum 
concentrations of alkaline phosphatase (ALP), aspartate 
aminotransferase (AST), blood urea nitrogen (BUN), and 


gamma glutamyltransferase (GGT) were elevated, and by days 
10-11 all animals had died. 


Tai Forest virus (TAFV) 

TAFV may cause severe disease and death in NHPs, as 
evidenced by the severe infection of one person after coming 
into contact with a dead chimpanzee, although the human 
patient ultimately survived (Guenno et al., 1995). To date, there 
have been no other reported infections, and thus this virus has 
not been the subject of extensive studies and there is a lack of 
animal models. 


SUMMARY 


Animal models are important for the study of filovirus 
pathogenesis, with the goal of testing the effectiveness of 
specific antivirals and vaccines. A variety of animal species are 
susceptible to WT and host-adapted filoviruses, and are thus 
used to test the efficacy of anti-filovirus products during non- 
outbreak situations. The Food and Drug Administration's "two- 
animal rule" states that the efficacy of a candidate product 
should be tested in “more than one animal species expected to 
react with a response predictive for humans" (FDA, 2017). 
Therefore, any vaccine or antiviral must demonstrate safety and 
efficacy in at least two animal models (a smaller animal model 
in addition to a NHP model) to have any hope of being licensed. 


The most important aspects of filovirus infection are 
summarized in Table 3. Parameters such as lymphopenia, liver 
damage, thrombocytopenia, coagulopathy, cytokine storm, 
rash, and hemorrhage are the chief hallmarks of filoviral 
diseases. Due to their advantages in terms of handling, cost, 
and access, mice and guinea pigs often serve as primary and 
secondary animal models for preliminary experiments and 
diverse genetic studies (such as knockout mice); however, as 
viruses can adapt to hosts and each host cannot recapitulate all 
hallmarks of filovirus diseases, mice and guinea pigs have 
limited translatability and predictability to human filoviral disease 
and thus are mostly used to screen potentially efficacious 
compounds. Currently, NHPs are the only animal species 
available that can recapitulate all stated aspects and thus are 
considered as the gold standard. Ferrets are the newest animal 
model to gain attention in filovirus research (Cross et al., 2016; 
Kozak et al., 2016; Kroeker et al., 2017b). In addition to being 
small, ferrets can accurately recapitulate human  filoviral 
diseases, including coagulopathy. Syrian golden hamsters can 
also recapitulate disease that closely mimics that in humans 
and NHPs, including coagulopathy. However, this model suffers 
from the need for host-adapted viruses and the lack of 
widespread reagents for studying immune responses. 
Therefore, ferrets and NHPs are likely the two best animal 
models at present for the testing of anti-filovirus compounds. 


Table 3 Clinical manifestations in different animal models of filovirus infections 





Immuno-competent mouse Immuno-comprimised mouse Guinea pig Syrian hamster Ferret NHP Human 
Host-adapted virus Yes Yes Yes Yes No* No* No* 
Lymphopenia Yes Yes Yes Yes Yes Yes Yes 
Liver damage Yes Yes Yes Yes Yes Yes Yes 
Thrombocytopenia Yes Yes Yes Yes Yes Yes Yes 
Coagulopathy No Unknown Yes Yes Yes Yes Yes 
Cytokine storm Yes Yes Unknown Yes Unknown Yes Yes 
Rash No No No Yes Yes Yes Yes 
Hemorrhage signs No Yes Unknown Yes Yes Yes Yes 


NHP: Non-human primates; No*: Wild-type nonadapted viruses. 


While this review has covered the different animal models 
available for each filovirus (Table 1), it is important to note that 
models against some viruses are more developed than others. 
There is an abundance of EBOV animal models, but very 
limited models available against TAFV. In addition, only ferrets 
have been used as a small animal model for BDBV. Immediate 
future goals should be the development and characterization of 
rodent animal models for the lesser known filoviruses. With 
experience gained from previous studies and well-established 
methods, the goal of establishing the necessary animal models 
for each known filovirus is attainable. 
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ABSTRACT 


Even with an effective vaccine, an estimated 240 
million people are chronically infected with hepatitis B 
virus (HBV) worldwide. Current antiviral therapies, 
including interferon and nucleot(s)ide analogues, 
rarely cure chronic hepatitis B. Animal models are 
very crucial for understanding the pathogenesis of 
chronic hepatitis B and developing new therapeutic 
drugs or strategies. HBV can only infect humans and 
chimpanzees, with the use of chimpanzees in HBV 
research strongly restricted. Thus, most advances 
in HBV research have been gained using mouse 
models with HBV replication or infection or models 
with HBV-related hepadnaviral infection. This review 
summarizes the animal models currently available for 
the study of HBV infection. 


Keywords: Hepatitis B virus; Animal model; Duck 
hepatitis B virus; Woodchuck hepatitis virus 


INTRODUCTION 


One-third of the world’s population have experienced hepatitis 
B virus (HBV) infection, with 240 million people chronically 
infected. Although hepatitis B vaccination programs have 
greatly reduced HBV infection rates in adolescents in China, 
there are still 93 million chronic HBV carriers (Chinese Society 
of Hepatology, Chinese Medical Association & Chinese Society 


of Infectious Diseases, Chinese Medical Association, 2015). 


Existing drug interferon and nucleoside (acid) analogues 
have low response rates and cannot effectively remove HBV; 
therefore, exploring new drugs or treatment programs is still an 
urgent problem in the field of HBV research. 

Animal models are promising tools for solving the above 
issue; however, such models are currently hampered 
because of the extremely narrow host range. HBV is 
the prototypic member of the Hepadnaviridae family, which 
contains closely-related hepatotropic, enveloped DNA viruses, 
including HBV, duck hepatitis B virus (DHBV), and woodchuck 
hepatitis virus (WHV). Chimpanzees are the only non-human 
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primate fully susceptible to HBV infection, though their use in 
HBV research is strongly restricted (Wieland, 2015). Several 
non-primate species can be infected with HBV, including 
the tree shrew (Tupaia belangeri), which also exhibits a 
close phylogenetic relationship with primates (Xiao et al., 
2017). In addition, primary Tupaia hepatocytes are suitable 
for HBV study, with sodium taurocholate co-transporting 
polypeptide (NTCP) recently identified using this system (Yan 
et al., 2012). Animals with non-HBV hepadnaviral infection 
have also been used in HBV research for many decades. 
Though DHBV is distantly related to HBV, primary duck 
hepatocytes and ducklings are easily available and have 
made great contribution in elucidating the unique replication 
mechanism of hepadnaviruses and in evaluating antiviral 
drugs (Mason, 2015). Furthermore, WHV infection in 
newborn American woodchucks (Marmota monax) can led 
to persistent infection and hepatocellular carcinoma (HCC) 
(Menne & Cote, 2007). Recently, we found the Chinese 
woodchuck (Marmota himalayana) to be highly susceptible to 
WHV infection (Wang et al., 2011). Woodchucks are widely 
used in exploring the pathogenesis of chronic hepadnavirus 
infection and HCC development, as well as in assessing 
antiviral drugs and immunotherapeutic approaches. Although 
the mouse is the best characterized and most convenient 
laboratory animal, it cannot be infected with HBV. However, 
after considerable effort over many decades, several mouse 
models with HBV replication and chimeric mouse models 
harboring human hepatocytes with HBV infection have been 
successfully established. In this review, we summarize the 
animal models currently available for the study of HBV infection. 
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ANIMAL MODELS WITH HBV INFECTION 


Chimpanzee model 

The chimpanzee is the only immunocompetent animal fully 
susceptible to HBV infection. After inoculation with serum from 
HBV patients, chimpanzees can develop acute and chronic 
HBV infection, accompanied by similar liver inflammation and 
cellular immune responses as patients with HBV infection 
(Wieland, 2015), which closely mimics the natural history of 
HBV infection in humans (Table 1). 

Chimpanzee models are essential for evaluating the safety 
and efficacy of HBV vaccines (McAuliffe et al., 1980), and have 
also been used to assess HBV therapeutic strategies (Eren 
et al., 2000). Recent studies have found that GS-9620, an 
oral agonist of toll-like receptor-7, can effectively reduce HBV 
titers in chimpanzees by inducing NK cell and virus-specific 
T cell responses (Lanford et al., 2013). Chimpanzees have 


also been used for exploring the pathogenesis of HBV infection. 


Using this model, a strong and polyclonal CD8* T cell response 
against HBV was determined to be a key factor in viral 
clearance through cytopathic and non-cytopathic pathways, 
with CD4* T cell priming crucial in the early phase of HBV 
infection (Asabe et al., 2009; Thimme et al., 2003). 

However, due to the ethical and availability limitations, 


chimpanzees are not routinely used in HBV research. 


Continuous efforts have been made to establish HBV infection 
in smaller non-human primates. Recently, cynomolgus 
monkeys from Mauritius were found with naturally occurring 
transmissible chronic HBV infection (Dupinay et al., 2013), 
and rhesus macaques can establish HBV infection after viral 
vector-mediated in vivo expression of human NTCP in its 
hepatocytes (Burwitz et al., 2017). 


Tupaia model 

The tree shrew (Tupaia belangeri) belongs to Tupaiidae, 
Scandentia, Mammalia, and displays a close phylogenetic 
relationship with primates. Recently, the whole genome 
sequence of the Chinese tree shrew implied that the nervous, 
immune, and metabolic systems were similar to those of 
humans (Xiao et al., 2017). To date, the tree shrew is the only 


non-primate animal found to be susceptible to HBV infection. 


Primary Tupaia hepatocytes are easily available compared with 
primary human hepatocytes, and therefore have been used 
in HBV research for many years (Table 1). NTCP, which is 
the cellular receptor responsible for HBV entry, was recently 
identified in the Tupaia model (Yan et al., 2012). Although 
HBV infection in neonatal tree shrews can lead to chronicity 
and pathological changes, including fibrosis, in vivo infection 
efficiency in Tupaia needs improvement (Wang et al., 2012; 
Yang et al., 2015). Recently, genotype A2 HBV isolates in 
Tupaia belangeri from Japan were found with higher HBV 
replication and chronicity rates, with the interferon response 
found to be impaired by HBV infection (Kayesh et al., 2017). 


Human chimeric mice 
The first human liver chimeric mouse model was developed 
in immunodeficient (Rag2-/-, SCID, SCID/beige) mice with the 
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urokinase-type plasminogen activator (uPA) transgene. The 
expression of the uPA gene can induce necrosis of hepatocytes, 
leading to subacute liver failure in young mice, and making it 
possible to transplant human hepatocytes into mouse livers. 
Transplantation of human hepatocytes into uPA-SCID mice 
results in a liver-humanized model with high human hepatocyte 
reconstitution rate and supportive of HBV and HCV infection 
(Dandri et al., 2001; Tsuge et al., 2005) (Table 1). 

A chimeric mouse model was constructed using FRG 
(Fah^/Rag2" /Il2rg^) mice, in which liver failure can be 
induced at will (Azuma et al., 2007). Results showed 
that fumarylacetoacetate hydrolase (Fah) deficiency led to 
liver failure, and double knock out of Rag2 and the 
gamma chain of the IL-2 receptor resulted in sufficient 
human liver cell repopulation in mice. Administration of 
2-(2-nitro-4-trifluoro-methylbenzoyl)-1,3-cyclohexanedione (NTBC) 
protected FRG mice from liver injury until NTBC withdrawal. 
Subsequent improvements, including refining transplantation 
procedures and increasing numbers of transplanted cells, 
increased liver cell repopulation rates to 95% in FRG mice 
(Bissig et al., 2010). Human liver cells repopulated in FRG mice 
support HBV and hepatitis C virus (HCV) infection (Table 1). 

The chimeric mouse model can be used for investigations 
on the biology of HBV infections, not only the distinct 
HBV genotypes but also naturally occurring variants or 
drug-resistant mutants. This model has also been used for 
studying interactions of HBV with the host and in testing novel 
antivirals. A new monoclonal antibody targeting a unique 
HBV surface antigen (HBsAg) epitope has been found to 
induce prolonged suppression of HBV in human hepatocytes 
transplanted into FRG mice (Zhang et al., 2016). Furthermore, 
NVR3-778, a novel capsid assembly modulator, was found 
to have high antiviral activity by reducing serum HBV DNA 
and HBV RNA levels in uPA/SCID mice with humanized 
livers (Klumpp et al., 2017). Two ribonuclease H inhibitors, 
a-hydroxytropolone and N-hydroxypyridinedione, have been 
found to suppress HBV replication in FRG human liver chimeric 
mice (Long et al., 2018). 

Because chimeric animals are genetically immunodeficient, 
they are unsuitable for the study of adaptive immunity or 
immunotherapy strategies. Researchers have used a variety 
of strategies to develop chimeric mice dually reconstituted 
with both immune cells and hepatocytes of human origin, 
including AFC8 (Washburn et al., 2011) and A2/NSG (Bility 
et al., 2014) mouse models. The former carries combined 
Rag2-/-/ll2rg-/- immunodeficiency genes, whereas the latter 
consists of NOD-SCID Il2rg-/- immunodeficient mice harboring 
the human HLA-A2 gene. Liver injury in both mouse models 
can be induced with inducible caspase-8 expression and 
CD95-activating antibodies, respectively. Liver cells and 
immune precursor cells from human fetal liver tissue have been 
transplanted into these two mouse models, with the AFC8 
model used to study HCV infection (Washburn et al., 2011) and 
the A2/NSG model shown to support persistent HBV infection, 
accompanied by liver inflammation and fibrosis (Bility et al., 
2014) (Table 1). 


Table 1 Characteristics of animal models for HBV study 





Virological characteristics 


Animal characteristics 











virus Vital Viral infection’ — cccDNA — Immunocompetent/ |. Class/Order/ References 
entry replication formation immunodeficient Family/Genus 
Mammalia/Primate/ 
Chimpanzee HBV Yes Infection Yes Immunocompetent No vt J ! Wieland, 2015 
Hominidae/Pan 
Xiao et al., 2017; 
Mammalia/Scandentia/ Yan et al., 2012; 
Tupaia HBV Yes Infection Yes Immunocompetent No K : 
Tupaiidae/ Tupaia Yang et al., 2015; 
Wang et al., 2012 
Bissig et al., 2010; 
Human Mammalia/Rodentia/ 
uma HBV Yes Infection Yes Immunodeficient Yes i i Dandri et al., 2001; 
chimeric mouse Muridae/Mus 
Tsuge et al., 2005 
Aves/Anseriformes/ Mala ge 
Duck DHBV Yes Infection Yes Immunocompetent No . Mason, 2015; 
Anatidae 
Schultz et al., 2004 
Cote et al., 2000; 
Mammalia/Rodentia/ Menne & Cote, 2007; 
Woodchuck WHV Yes Infection Yes Immunocompetent No MS 
Sciuridae/Marmota Roggendorf et al., 2007; 
Wang et al., 2011 
HBV i Immunocom nt, Mammalia/Rodentia/ VEM 
transgenic HBV No Replication No mmunocomperent Yes ; i Guidotti et al., 1995 
mouse model tolerant to HBV Muridae/Mus 
Huang et al., 2006, 2012; 
HBV transfected Mammalia/Rodentia/ 
HBV No Replication No Immunocompetent Yes : Peng et al., 2015; 
mouse model Muridae/Mus 


ANIMAL MODEL WITH NON-HBV HEPADNAVIRAL 
INFECTION 


Duck model 

Mason et al. (1980) found a DNA virus similar to human 
HBV in the serum of ducks, named DHBV (Table 1) and 
subsequently classified as Avihepadnavirus. Similar to HBV 
infection in humans, age and inoculum dose can influence the 
outcome of DHBV infection. For example, 7-day-old ducklings 
can exhibit persistent infection when inoculated with high-titer 
DHBV inoculum (Cote et al., 2000). The rates of natural 
HBV infection can also vary in distinct duck species, ranging 
from 8.75% to 17.80%, and the rates of experimental DHBV 
infection can differ in ducklings inoculated via intraperitoneal 
or intravenous injection (6596 vs. 
Duck models with persistent DHBV infection are widely used 
to evaluate anti-HBV drugs and therapeutic strategies; for 
example nucleot(s)ide analogs commonly used in clinicals 
such as entacavir (ETV) (Marion et al., 2002), nucleocapsid 
assembly inhibitors (Campagna et al., 2013), nucleic acid 
polymer REP 2139 (Quinet et al., 2017), and combined 
immunotherapy (Le Guerhier et al., 2003). This animal model 


has also been used to study the pathogenesis of HBV infection. 


A lack of innate immunity rather than adaptive immunity in 
early phase (within 5 days after DHBV inoculation DHBV 
infection is important in DHBV persistence in newborn ducks 
(Tohidi-Esfahani et al., 2010). Residual cccDNA pools are 
detectable and unchanged in ducks following resolution of 
acute HBV infection, even after ETV treatment (Reaiche et 


8096) (Hao et al., 2012). 


Yang et al., 2002 


al., 2010). After continuous DNA vaccination combined with 
interleukin (IL)-2 and interferon (IFN)-y plasmid injections, only 
trace cccDNA is detectable in duck liver, whereas inoculation 
of naive ducks with these liver homogenates can cause high 
levels of DHBV infection (Saade et al., 2013). 


Woodchuck model 

In 1978, WHV was found in a colony of American woodchucks 
in Penrose Zoo, Philadelphia, with high rates of chronic 
hepatitis and HCC observed (Summers et al., 1978), and was 
subsequently identified as Orthohepadnavirus. WHV is very 
similar to HBV not only in terms of virological characteristics, 
but also the natural history of infection and immune response 
against viral infection (Cote et al., 2000) (Table 1). Therefore, 
it is widely used to study the pathogenesis of HBV infection 
(Fletcher et al., 2013; Menne et al., 1998), and to evaluate 
antiviral drugs, prophylactic vaccines, and immunotherapeutic 
strategies (Colonno et al., 2001; Kosinska et al., 2013; 
Roggendorf et al., 2007). 

We found the Chinese woodchuck (Marmota himalayana), a 
relative of the American woodchuck, to be highly susceptible 
to WHV, suggesting potential as a new animal model for 
HBV infection (Wang et al., 2011). However, the lack of 
genomic information on the Chinese woodchuck has strongly 
limited its application in HBV research. Recently, however, 
the transcriptome sequences of this species have been fully 
sequenced and analyzed, with high sequence similarity found 
between the Chinese and American woodchucks (Liu et al., 
2016). To improve the quality of the Chinese woodchuck 
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model, we inoculated animals with different dosages of WHV 
inoculum (104, 106, or 108 genome equivalents of WHV7), 
and found that different doses had a strong influence on the 
course of WHV infection (Wang et al., 2014); we also inoculated 
animals from different counties (i.e., Tongren, Tongde, Guide, 
and Wulan), and found varied susceptibilities to WHV infection 
(unpublished data). For better application of the Chinese 
woodchuck model, breeding colonies have been established in 
Qinghai and Xinjiang in northwest China. Furthermore, over 
20 genes from woodchucks have been analyzed, including 
IL-6 receptor, IL-12 receptor, IL-18 and its receptor, IL-21 
receptor, tumor growth factor-Bs and their receptors, CD40, 
CD160, CD244, Tim-3, galectin-9, retinoic acid-inducible gene 
l, melanoma differentiation-associated gene 5, AIM-2, IFI16, 
RelA, and interferon stimulated gene (Fan et al., 2012; Jiang 
et al., 2012; Lu et al., 2008; Wang et al., 2016; Yan et al., 
2016; Yang et al., 2013). These efforts have improved the 
Chinese woodchuck model for application in HBV research. An 
alternative strategy for HBV post-exposure prophylaxis based 
on nucleotide analogues (NA) was evaluated using this model, 
with NA alone or in combination with a DNA vaccine found to 
completely prevent viremia after WHV inoculation, but induce 
partial or complete protective immunity, respectively (Wang et 
al., 2014). 


MOUSE MODEL SUPPORTING HBV REPLICATION 


HBV transgenic mouse model 

Since the 1980s, HBV transgenic mice selectively expressing 
HBV proteins, including HBsAg, e antigen, or x protein, 
have been used to study the role of these proteins. In 
1995, a transgenic mouse model for the full HBV genome 
was established, which produced infectious HBV virions 
morphologically indistinguishable from human-derived virions 
(Guidotti et al, 1995) (Table 1). Although HBV is 
immune-tolerant and cannot be spontaneously cleared as its 
genome is integrated in the mouse genome, this model has 
made important contributions in elucidating the pathogenesis 
of HBV infection: for example, no liver damage has been 
observed in HBV transgenic mice, suggesting that HBV itself 
is not cytopathic; and the adoptive transfer of HBV-specific 
cytotoxic T cells can lead to liver damage in HBV transgenic 
mice and transient clearance of HBsAg through lysis of 
HBsAg-expressing cells, whereas viral replication can also 
be suppressed by non-cytolytic mechanisms involving IFN-o 
and tumor necrosis factor-a (TNF-o) (Guidotti et al., 1996; 
Heise et al., 1999). HBV transgenic mice have also been 
used for testing antiviral drugs, such as lamivudine, adefovir 
dipivoxil, ETV, and HBV-targeting siRNAs (Ebert et al., 2011; 
McCaffrey et al., 2003). In China, HBV transgenic mice 
have also been established and applied in HBV research, 
including: (1) evaluating a DNA vaccine by fusion of the 
granulocyte-macrophage colony stimulating factor gene to the 
HBV-S gene (Qing et al., 2010); (2) analyzing hepatocyte 
proteomics and seroproteomics in HBV-transgenic mice (Ding 
et al., 2009); and (3) exploring the role of NK cells in CCL4 
accelerated liver fibrosis (Jin et al., 2011). 
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HBV transfected mouse model 

Viral vectors based on adenoviruses or adeno-associated viruses 
containing HBV genomes efficiently transduce hepatocytes in 
immunocompetent mice. After intravenous injection of high 
doses of HBV containing adenoviral vectors (Ad-HBV), HBV 
replication lasts for only three months because transduction of 
ad-HBV often results in a strong immune response against the 
adenovirus itself. This can be overcome by the injection of low 
doses of Ad-HBV (Huang et al., 2012) (Table 1). 

Injecting a large amount of liquid containing naked DNA 
into mice by the tail vein over a short period of time terms 
as hydrodynamic injection, which can effectively transfer 
exogenous genes into liver cells. Hydrodynamic injection 
was found to successfully deliver replication-competent HBV 
genomes into the livers of immunocompetent mice, with viremia 
peaking after 6 d and rapidly declining thereafter (Yang et 
al., 2002). It should be noted that plasmid backbone, mouse 
strain, and sex can affect HBV replication duration in mice. 
HBV replication was found to last for up to six months using 
the pAAV/HBV1.2 plasmid and C57BL/6 mice (Huang et al., 
2006), and up to 46 weeks using the pAAV/HBV1.2 plasmid 
and C3H/HeN mice (Peng et al., 2015) (Table 1). This animal 
model can be used to analyze immune responses during acute 
or persistent HBV replication, clarify the biology of different 
HBV genotypes, variants, or mutants in vivo, and evaluate 
antiviral compounds. Using this model, we discovered that Poly 
(I: C) cleared HBV through IFN-dependent pathways (Wu et al., 
2014). Recently, one genotype B HBV isolate, designated BPS, 
persisted for 38 weeks in ~50% of mice after hydrodynamic 
injection, with IL-21 and IL-33 found to be potent regulators of 
HBV clearance (Shen et al., 2017). 

HBV cccDNA plays a major role in HBV persistence and 
is responsible for the relapse of viral activity after antiviral 
treatment with polymerase inhibitors in chronically infected 
patients, while is not built in not only HBV transgenic mice but 
also HBV transfected mice. To overcome this large obstacle 
and closely mimic chronic HBV infection in humans, Li et 
al. (2018) established a surrogate mouse model with HBV 
CCCDNA persistence; using a replication-defective recombinant 
adenoviral vector, recombinant HBV cccDNA was successfully 
delivered into the mouse liver and was intrinsically stable, 
resulting in HBV persistence throughout the experiment (262 
weeks) in transgenic mice expressing Cre recombinase under 
the albumin promoter. 


SUMMARY 


The above animal models with HBV replication or hepadnavirus 
infection provide more options for HBV research. The 
identification of NTCP as an HBV receptor is a milestone in 
the field of HBV research. The NTCP-transfected hepatoma 
cell line is the first human cell line supporting HBV infection, 
while NTCP transgenic mice support HBV infection was 
ongoing. Transfected, transduced, or transgenic mice with 
HBV genomes can only support HBV replication, with the 
HBV life cycle incomplete due to the lack of viral entry, 
cccDNA formation, and viral spread. As human liver cells 


supporting HBV infection are transplanted into immunodeficient 
mice, these mice exhibit obvious immunodeficiency and their 
maintenance systems are very complex. The chimpanzee is 
an optimal model for HBV infection, but is strongly restricted 
and not easily available for HBV research. Tree shrews can be 
infected with HBV, but the in vivo system requires considerable 
improvement. For the duck and woodchuck models, the 
differences between HBV and other hepadnaviruses should 
be considered. In addition, except for mice, the above 
models are not inbred and detection reagents are not readily 
available. Thus, all current animal models have specific 
limitations. Therefore, researchers need to carefully interpret 
their results from animal studies, and validation of their findings 
in multiple systems should be encouraged. 
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ABSTRACT 


The Angolan strain of Marburg virus (MARV/Ang) 
can cause lethal disease in humans with a case 
fatality rate of up to 90%, but infection of 
immunocompetent rodents do not result in any 
observable symptoms. Our previous work includes 
the development and characterization of a 
MARV/Ang variant that can cause lethal disease in 
mice (MARV/Ang-MA), with the aim of using this tool 
to screen for promising prophylactic and therapeutic 
candidates. An intermediate animal model is needed 
to confirm any findings from mice studies before 
testing in the gold-standard non-human primate 
(NHP) model. In this study, we serially passaged the 
clinical isolate of MARV/Ang in the livers and 
spleens of guinea pigs until a variant emerged that 
causes 100% lethality in guinea pigs (MARV/Ang- 
GA). Animals infected with MARV/Ang-GA showed 
signs of filovirus infection including lymphocytopenia, 
thrombocytopenia, and high viremia leading to 
spread to major organs, including the liver, spleen, 
lungs, and kidneys. The MARV/Ang-GA guinea pigs 
died between 7-9 days after infection, and the LDso 
was calculated to be 1.1x10-' TCIDso (median 
tissue culture infective dose). Mutations in 
MARV/Ang-GA were identified and compared to 
sequences of known rodent-adapted MARV/Ang 
variants, which may benefit future studies 
characterizing important host adaptation sites in the 
MARV/Ang viral genome. 


Keywords: Marburg virus; Guinea pig; Animal model; 
Pathogenesis; Host adaptation 
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INTRODUCTION 


Marburg virus (MARV) is a pathogen belonging to the family 
Filoviridae, genus Marburgvirus, and is classified as a Risk 
Group 4 Pathogen (Phac-Aspc.Gc.Ca, 2011) MARV is 
considered to be extremely pathogenic in humans. Initial 
symptoms of infection are nonspecific with high grade fever, 
chills, headache and muscle aches followed by nausea, 
vomiting, diarrhea and abdominal pain. Progression of disease 
to advanced stages is characterized by the appearance of 
maculopapular rash, increased vascular permeability, 
neurological signs, severe dehydration, diffuse coagulopathy, 
hemorrhage, shock and multi-organ failure leading to death 
within 7—10 days after symptoms appear (Mehedi et al., 2011). 
The few patients that survive the infection enter a prolonged 
convalescence phase but are beset with complications 
including ocular symptoms, joint pain, orchitis and psychosis, 
among others (Mehedi et al., 2011). Additionally, MARV is 
known to persist in some survivors and can be spread via 
semen (Brainard et al., 2016), potentially resulting in secondary 
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cases of infection. There are currently no specific preventatives 
or cures approved for treating MARV infections. 

The first MARV outbreak was reported during 1967 in 
Marburg, Germany, in which 7 deaths from 31 total infections 
were recorded. The source of the virus was eventually traced to 
infected monkeys imported from Uganda (Kissling et al., 1968; 
Smith et al., 1967). A total of 12 outbreaks have been recorded 
so far, with the last occurring in Uganda during 2014 (CDC Gov, 
2014). The largest and deadliest MARV outbreak occurred 
during 2004-2005 in Angola when 227 fatal cases out of 252 
infections were reported, for a case fatality rate (CFR) of ~90% 
(Towner et al., 2006). Experimental infection of cynomolgus 
macaques with the Angolan strain of MARV (MARV/Ang) 
suggests that this virus appears to cause a more rapid onset 
and severe outcome of infection, compared to other MARV 
strains (Fernando et al., 2015). Live MARV has been isolated 
from bats (Towner et al., 2009) and several outbreaks have 
been associated through contact with wild bats inside mines 
and caves (Adjemian et al., 2011; CDC, 2009; Timen et al., 
2009), strongly suggesting that bats play a major role in 
transmission of the virus to humans. 

Due to their low unit costs and ease of handling, rodents are 
particularly desirable for initial studies in vivo, especially for the 
screening of candidate antivirals against infectious diseases. 
However, immunocompetent rodents are not susceptible to 
infection with wild-type MARV isolates. To address this issue, we 
had previously developed a variant of MARV/Ang via serial 
passaging in the livers of severe combined immunodeficiency 
(SCID) mice (Qiu et al., 2014). This novel, mouse-adapted variant 
(MARV/Ang-MA) caused uniformly lethal disease when injected 
intraperitoneally (i.p.) into BALB/c and C56BL/6 mice, replicating to 
high levels in the blood, intestines, kidneys, lungs, brain, spleen, and 
livers. Other findings include lymphocytopenia, thrombocytopenia 
and marked damage to the liver and spleen (Qiu et al., 2014). The 
MARV/Ang-MA model is currently being used to assess the 
potential of candidate drugs against MARV (Warfield et al., 2017). 

Non-human primates (NHPs) are widely considered to be the 
best animal model available for recapitulating human MARV 
disease, and any candidate antivirals must show efficacy in 
NHPs before clinical trials. However, it is prudent to develop 
and characterize an alternate small animal model to confirm 
any findings from mice before testing in NHPs, as these 
animals are costly, laborious to handle, and require specialized 
housing. To address this need, we developed a MARV/Ang 
variant adapted to guinea pigs (MARV/Ang-GA) by sequential 
passaging in the livers and spleens of these animals. 
MARV/Ang-GA was then characterized with regards to median 
lethal dose (LDso), viremia/biodistribution at various times after 
infection. The full length viral genome of MARV/Ang-GA was 
also sequenced to identify any mutations arising from the 
passaging process, and compared to the wild-type MARV/Ang 
as well as available sequences of rodent-adapted MARV/Ang to 
elucidate potential “hotspots” for adaptation mutations. 


MATERIALS AND METHODS 


Cells and viruses 
The wild-type, clinical virus isolate Marburg virus/H.sapiens- 


tc/AGO/2005/Angola, or MARV/Ang, was isolated from a 
patient during the 2005 Marburg outbreak in Uige, Angola. The 
adapted virus is called Marburg  virus/NML/C.porcellus- 
lab/AGO/2005/Ang-GA-P2, or MARV/Ang-GA. Stocks of 
MARV/Ang and MARV/Ang-GA were grown in T-150 flasks 
(Corning, USA) of Vero E6 cells (ATCC, USA) for these studies. 


Animals and virus passaging 

Outbred female, 4—6 week old strain Hartley guinea pigs 
between 200-250 g in weight were used for passaging 
experiments. Two guinea pigs were injected i.p. with ~10° 
plaque forming units (pfu) of MARV/Ang in 1 mL of Dulbecco's 
modified Eagle's medium (DMEM), supplemented with 296 
heat-inactivated fetal bovine serum (FBS) (Sigma, Canada). 
Serial passaging of the virus was then performed as follows. 
Livers and spleens were removed from euthanized animals at 7 
days post infection (dpi), pooled, and homogenized by grinding 
the livers against a steel mesh with a plastic plunger. 
Splenocytes and hepatocytes were suspended in 10 mL of 
phosphate-buffered saline (PBS) and centrifuged at 400 r/min 
for 5 min. The supernatant was passed through a 70 um cell 
strainer, and the cell pellet was mechanically homogenized in 1 
mL DMEM with sterile steel beads in a tissue homogenizer. 
After centrifugation at 400 r/min for 5 min, the supernatant from 
this cell pellet was also passed through the 70 um cell strainer. 
Two naive guinea pigs were injected i.p. with 1 mL of the 
homogenate, and the process was repeated again at 7 dpi until 
MARV/Ang-GA was produced after 9 passages. MARV/Ang-GA 
was not plaque purified. 

All animal work was approved by the Animal Care Committee 
(ACC) at the Canadian Science Center for Human and Animal 
Health (CSCHAH), and performed in accordance with the 
guidelines from the Canadian Council on Animal Care (CCAC). 
Studies involving live, infectious virus were performed under 
biosafety level 4 (BSL-4) conditions at the National Microbiology 
Laboratory located in Winnipeg, Canada. 


LDso determination 

Serial 10-fold dilutions of MARV/Ang-GA from 2.32x10-2— 
2.32x10? TClDso (median tissue culture infective dose) was 
administered per guinea pig (n=3 per group). Animals were 
weighed daily and monitored for survival and weight loss. 


Serial sampling of infected guinea pigs 

Guinea pigs were infected with 1 000xLDs of MARV/Ang-GA, 
or an equivalent dose of MARV/Ang. Four animals from each 
group were euthanized at 0, 3, 5 and 7 dpi for necropsy. Blood 
was sampled from anesthetized guinea pigs by cardiac puncture. 
Whole blood was collected in K2 EDTA plus blood collection 
tubes (BD Biosciences, Canada) for blood counts and stored in 
lysis buffer AVL (Qiagen, USA) for determination of viremia, and 
lithium heparin blood collection tubes (BD Biosciences, Canada) 
for blood biochemistry. Major organs including livers, spleen, 
kidneys and lungs were harvested from these animals at the 
same timepoints and stored in RNAlater for subsequent analysis. 


Blood counts and blood biochemistry 
Complete blood counts and blood biochemistry analyses 
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were performed with a VetScan HM5 hematology system or 
a VetScan VS2 analyzer, respectively (Abaxis Veterinary 
Diagnostics, USA). 


Determination of viral titers by reverse transcription- 
quantitative PCR (RT-qPCR) 

RNA was extracted from whole blood using the QlAamp viral 
RNA minikit (Qiagen, USA) or from tissue homogenates with an 
RNeasy minikit (Qiagen, USA), following manufacturer 
instructions. The amount of viral genome was quantified via 
RT-qPCR on the ABI StepOnePlus, using the LightCycler 
480 RNA master hydrolysis probe kit (Roche, Canada) with 
the RNA polymerase as the target gene. Reaction conditions 
were as follows: 63 °C for 3 min, 95 °C for 30 s, and then 45 
cycles of 95 °C for 15 s and 60 °C for 30 s. Primer and 
probe sequences were as follows: forward, 5'-GCAAAAG 
CATTCCCTAGTAACATGA-3'; reverse, 5-CACCCCTCACTA 
TRGCGTTYTC-3’; probe, 6-carboxyfluorescein [FAM]-TGGC 
ACCAYAATTCAGCAAGCAT AGG-black hole 
quencher[BHQ]. 


Histopathology and immunohistochemistry 

Tissues were fixed in 10% neutral phosphate buffered formalin, 
routinely processed, sectioned at 5 um and stained with 
hematoxylin and eosin (HE) for histopathologic examination. 
For immunohistochemistry (IHC), paraffin tissue sections were 
quenched for 10 min in aqueous 3% hydrogen peroxide. 
Epitopes were retrieved using an in-house glyca retrieval 
solution, in a Biocare Medical Decloaking Chamber. The 
primary antibody applied to the sections was an in-house anti- 
MARV mouse monoclonal antibody (3H1). It was used at a 
1:4 000 dilution for 30 min. They were then visualized using a 
horse radish peroxidase labelled polymer, Envision®+system 
(anti-mouse) (Dako, USA) and reacted with the chromogen 
diaminobenzidine (DAB). The sections were then counter 
stained with Gill’s hematoxylin. 


Statistical analysis 

Statistical analysis was performed with the two-way ANOVA 
with Bonferroni post-tests. Statistical significance was set at 
*: P<0.05, **: P<0.01, or ***: P<0.001. 


Sequencing of MARV/Ang-GA 

RNA was extracted from MARV/Ang and MARV/Ang-GA stocks 
using a Qiagen QlAmp viral RNA minikit (Qiagen, USA), 
following manufacturer instructions. Sanger sequencing of 
MARV/Ang-GA is carried out using the same protocol as 
described in a previous publication (Qiu et al., 2014). 


Nucleotide sequence accession number 

The sequence of MARV/Ang-GA was deposited in GenBank 
(Accession No. MF939097). 

RESULTS 


LDso determination of MARV/Ang-GA in guinea pigs 
Groups of three guinea pigs were inoculated i.p. with serial 10- 
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fold dilutions of MARV/Ang-GA, ranging from 2.32x10—— 
2.32x10? TCIDso, in order to assess the LDso of this adapted 
virus. Survival and weight loss was monitored for 15 days after 
challenge. The animals that were challenged with 2.32x109— 
2.32x10? TCIDso of MARV/Ang-GA all succumbed to disease on 
7-8 dpi (the mean time to death for the 2.32x10? TCIDso group 
was 7.3+0.6 dpi) with close to 20% weight loss (Figures 1A, B). 
The animals that were challenged with 2.32x10-' TCIDso of 
MARV/Ang-GA died on 8-9 dpi with almost 20% weight loss 
(Figures 1A, B) The animals given 2.32x10? TClDso of 
MARV/Ang-GA all survived the challenge with no weight loss 
observed (Figures 1A, B). Using the linear regression from the 
logarithmic value of the challenge titer and the fatality rate, the 
LDso was estimated to be approximately 1.1x10-' TCIDso. 
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Figure 1 LDso determination of MARV/Ang-GA in guinea pigs 
A-B: Survival (A) and weight loss (B) of guinea pigs infected 
intraperitoneally with serial 10-fold dilutions of MARV/Ang-GA ranging 
from 2.32x107—2.32x10? TCIDso per animal. 


Clinical overview of changes in guinea pigs after infection 
with MARV/Ang or MARV/Ang-GA 

Groups of three guinea pigs were then challenged i.p. with 
1 000xLDso (or 110 TCIDso) of MARV/Ang-GA, or an equivalent 
dose of MARV/Ang. As expected, weight loss and clinical 
symptoms were only observed in the MARV/Ang-GA animals, 
but not in the MARV/Ang group (Supplementary Figure 1A, B). 
A complete blood count as well as biochemical analysis was 
performed from whole blood drawn from infected animals at 0, 
3, 5 and 7 dpi. A drop in the counts of white blood cells and 
lymphocytes was observed in MARV/Ang-GA animals at 5 dpi 
(Figures 2A, B), whereas the decrease in lymphocyte 


percentage was apparent and statistically significant from 3—7 
dpi (Figure 2C). Additionally, a significant decrease in platelet 
count was observed in MARV/Ang-GA animals at 3—5 dpi, 
whereas this decrease was more gradual and less pronounced 
in MARV/Ang animals (Figure 2D). Advanced to terminal 
disease in MARV/Ang-GA guinea pigs can be observed by 
statistically significant increases in the following markers: 
albumin, alkaline phosphatase, alanine aminotransferase, total 
bilirubin, blood urea nitrogen, creatinine and globulin, while 
these markers stayed relatively constant in MARV/Ang animals 
(Figures 3A—G). Additionally, a decrease in total protein levels 
was observed in MARV/Ang-GA guinea pigs at 7 dpi, but not in 
MARV/Ang animals (Figure 3H). 


Viremia and biodistribution of MARV/Ang or MARV/Ang- 
GA in guinea pigs 

The amount of MARV/Ang or MARV/Ang-GA in the blood and 
selected major organs of infected guinea pigs were quantified 
by RT-qPCR at 0, 3, 5 and 7 dpi. Animals infected with 
MARV/Ang showed evidence of viral replication, with viremia of 
-10? rising to ~10* genome equivalents per mL of blood 
(GEQ/mL blood) from 3-7 dpi, whereas those infected with 
MARV/Ang-GA showed much higher levels of replication, with 
~10° rising to 10’ GEQ/mL blood from 3-7 dpi (Figure 4A). 
MARV/Ang was found at substantial levels in the livers and 
spleens of guinea pigs, with peak values of ~10* and ~10° GEQ 
per gram of tissue (GEQ/g tissue), respectively, whereas for 
MARV/Ang-GA the peak values in the liver and spleen were 
~10’ GEQ/g tissue (Figures 4B, C). While substantial levels of 
MARV/Ang-GA was found in the kidneys and lungs of infected 
animals with peak values of -109 GEQ/g tissue in both organs 
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at 7 dpi, MARV/Ang is not as readily detectable in the kidneys 
or lungs, with peak values of ~10° GEQ/g tissue for both organs 
(Figures 4D, E). 


Histopathology and immunohistochemistry 

In histopathology results, lesions were observed in all 
tissues examined from infected animals but no lesions were 
observed in uninfected control animals (Figures 5A-H). In 
the lungs the alveolar walls were expanded by infiltration of 
inflammatory cells and there was extensive individual cell 
necrosis. Throughout the liver there was degeneration of 
hepatocytes often with mineralization as well as portal 
inflammation and intracytoplasmic inclusion bodies. In the 
spleen there was massive necrosis of cells within both the 
red and white pulp areas. Mild lesions in the kidneys were 
characterized by scattered individual degeneration of cells 
within glomeruli. In  immunohistochemistry results viral 
antigen was detected in all infected animals including lung 
alveolar walls, sinusoids/ hepatocytes of the liver, red and 
white pulp of the spleen and within the renal glomeruli 
however no immunostaining was observed in uninfected 
control animals (Figures 6A—-H). 


Mutations in MARV/Ang-GA 

The adapted virus MARV/Ang-GA was then sequenced and 
aligned with the publicly available sequence of MARV/Ang 
(GenBank accession No. DQ447653), in order to identify any 
mutations that arose and became fixed through passaging. A 
total of eight mutations were identified: two mutations in non- 
coding regions (base pairs 2 931 and 18 713), two mutations in 
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Figure 2 Complete blood counts of guinea pigs challenged with 1 000xLDso of either MARV/Ang-GA or MARV/Ang 
A: White blood cells (WBC); B: Lymphocytes (LYM); C: Lymphocyte percentage (LYM%); D: Platelet (PLT); *: P<0.05; **: P<0.01; ***: P< 0.001. 
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Figure 3 Blood biochemistry of guinea pigs challenged with 1 000xLDso of either MARV/Ang-GA or MARV/Ang 
A: Albumin (ALB); B: Alkaline phosphatase (ALP); C: Alanine aminotransferase (ALT); D: Total bilirubin (TBIL); E: Blood urea nitrogen (BUN); F 
Creatinine (CRE); G: Globulin (GLOB); H: Total protein (TP); *: P«0.05; **: P«0.01; ***: P«0.001. 
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Figure 4 Viremia and viral biodistribution of guinea pigs challenged with 1 000xLDso of either MARV/Ang-GA or MARV/Ang 
The amount of viral RNA was quantified by RT-qPCR in the blood (A), liver (B), spleen (C), kidney (D) and lungs (E) of the infected animals; the dashed 
lines represent the limit of detection with the RT-qPCR assay; *: P<0.05; **: P<0.01; ***: P<0.001. 


the polymerase (L) coding region (base pairs 13115 and 
17 249) resulting in silent mutations at the amino acid level, one 
mutation in the viral protein 40 (VP40) coding region (base pair 
4 735) resulting in a N56K change at the amino acid level, and 
three mutations in the viral protein 24 (VP24) coding region (base 
pairs 10 402, 10 853 and 10 931) resulting in changes of V66l, 
L216S and N242S, respectively, at the amino acid level (Table 1). 


DISCUSSION 


Using the well-established method of passaging wild-type 
filovirus isolates in the livers and spleens of rodents, we 
successfully generated and characterized a guinea pig model 
against MARV/Ang, currently the most lethal strain of MARV 
to humans. Animals infected with MARV/Ang-GA developed 
severe disease characterized by weight loss, decrease in 
white blood cells and platelets, as well as increased 
levels/activities of markers for renal, hepatic and pancreatic 
functions, indicating multi-organ failure. High viremia and 
systemic spread of the virus was detected during the course 
of the disease. 


Another research group had also reported on the 
development of a guinea pig animal model based on 
MARV/Ang (Cross et al., 2015). In their model, animals infected 
with 5 000 pfu of the guinea pig-adapted virus experience fever 
and weight loss resulting in death between 8—10 dpi. Gross 
pathology, histopathological and immunohistochemical studies 
show massive hepatic and splenic damage. The hematology 
and biochemical findings show lymphocytopenia and 
thrombocytopenia, but marked increases in liver-associated 
enzyme and pro-inflammatory cytokine levels. Additionally, 
they also show that guinea pigs infected with the adapted 
virus have coagulopathy, as evidenced by increased 
prothrombin and activated partial thromboplastin time, as 
well as decreased thrombin times and protein C activity 
(Cross et al., 2015). 

Since the virus passaging by Cross et al. (2015) was 
performed completely independently from our study, this 
provided a rare opportunity to compare viral genome 
sequences between their guinea pig-adapted virus versus our 
MARV/Ang-GA. The purpose was to see if any common 
mutations exist between the two adapted viruses, which may 
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Figure 5 Histopathology findings in uninfected versus GA-MARV infected guinea pigs 

Lung (A), liver (C), spleen (E) and kidney (G) from uninfected control guinea pigs were normal and no significant lesions were detected. In 
infected guinea pigs lesions were observed in all tissues; B: In the lung, the alveolar walls contained numerous inflammatory cells and there 
were many degenerating/necrotic cells throughout (arrows); D: Liver lesions included hepatocyte degeneration/necrosis often with 
mineralization (thin arrows) and numerous intracytoplasmic viral inclusion bodies were observed (inset); F: In the spleen, there was massive 
necrosis in both the red pulp (*) and white pulp (arrow) areas; H: In the kidneys, there were degenerating cells with pyknotic nuclei scattered 
throughout the glomeruli (arrows). Scale bars: 20 um. 
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Figure 6 Immunohistochemistry findings in uninfected versus GA-MARV infected guinea pigs 


No positive immunostaining for MARV viral antigen was observed in lung (A), liver (C), spleen (E) or kidney (G) from uninfected control guinea pig. 
Abundant MARV antigen was detected in lung (B), liver (D), spleen (F) and kidney (H, arrows point to cells infected within glomeruli) from infected 


guinea pigs. Scale bars: 20 um. 
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give indications of the relative importance of each mutation. 
To our surprise, the two viruses were remarkably similar. All 
mutations in the coding and non-coding regions were 
identical, save for the N242S mutation in VP24 described in 
our study that was not present in their variant, and a mutation 
in a non-coding region (base pair 19 105) that was described 
in their study but undetected in our viral stock (Table 1). This 
suggested that the common mutations between our two viral 
Stocks were the critical changes needed to attain virulence in 


guinea pig; however, it should be noted that some of the 
mutations may have been introduced during the preparation of 
the virus stock in VeroE6 cells, as previously demonstrated 
with a mouse-adapted MARV/Ang (Wei et al., 2017). It will be 
interesting to sequence the virus from the last in vivo passage 
from the guinea pigs and compare with the sequence of our 
viral stock. Additionally, in-depth studies using the reverse 
genetics system for MARV will be needed to elucidate the 
impact of each mutation. 


Table 1 Summary of mutations in MARV/Ang-GA, compared to MARV/Ang 








icd Coding (gene) MARV/Ang MARV/Ang-GA (Cross et al., MARV/Ang-GA (this Amino acid 

position (bp) or non-coding 2015) study) change? 

2 931 Non-coding U A A N/A 

4 735 VP40 U A A N56K 

10 402 VP24 G A A V66I 

10 853 VP24 U C C L216S 

10 931 VP24 A A G N2428 

iem : K e T PAM 

17 249 L U A A Silent 
mutation 

18 713 Non-coding C A A N/A 

19 105 Non-coding A U A N/A 


N/A: not applicable. 


Another important finding is that although non-coding regions 
and silent mutations are not expected to theoretically have an 
impact on virus evolution, the similarity of the mutations in these 
locations between the two viruses, despite independent 
passaging, suggests that the non-coding regions and silent 
mutations may play a larger-than-expected role. It will be 
important to also sequence the viral homogenates during each 
passage to observe when the mutations described arise and 
become fixed, and whether the mutations always appear in the 
same order. These studies will help determine whether viral 
evolution and host adaptation is based on a set template plan, 
and may give important clues on predicting viral evolution and 
host adaptation during subsequent MARV outbreaks. 
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ABSTRACT 


Parasites can increase infection rates and 
pathogenicity in immunocompromised human 
immunodeficiency virus (HIV) patients. However, in 
vitro studies and epidemiological investigations 
also suggest that parasites might escape 
immunocompromised hosts during HIV infection. 
Due to the lack of direct evidence from animal 
experiments, the effects of parasitic infections on 
immunocompromised hosts remain unclear. Here, 
we detected 14 different parasites in six northern 
pig-tailed macaques (NPMs) before or at the 50th 
week of simian immunodeficiency virus (SIV) infection 
by ELISA. The NPMs all carried parasites before viral 
injection. At the 50th week after viral injection, the 
individuals with negative results in parasitic detection 
(i.e., 08247 and 08287) were characterized as the 
Parasites Exit (PE) group, with the other individuals 
(i.e., 09203, 09211, 10205, and 10225) characterized 
as the Parasites Remain (PR) group. Compared with 
the PR group, the NPMs in the PE group showed 
higher viral loads, lower CD4* T cells counts, and 
lower CD4/CD8 rates. Additionally, the PE group had 
higher immune activation and immune exhaustion of 
both CD4* and CD8* T cells. Pathological observation 
showed greater injury to the liver, cecum, colon, 


spleen, and mesenteric lymph nodes in the PE group. 


This study showed more seriously compromised 
immunity in the PE group, strongly indicating that 
parasites might exit an immunocompromised host. 
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INTRODUCTION 


Acquired Immune Deficiency Syndrome (AIDS), induced by 
human immunodeficiency virus (HIV), was first recognized in 
1981 and has emerged from initially being a concern among 
high-risk groups to representing a worldwide pandemic (Jeang 
et al., 2007). According to estimates by WHO and UNAIDS, 
36.7 million people were living with HIV globally at the end 
of 2015. That same year, approximately 2.1 million people 
became newly infected and 1.1 million died of HIV-related 
causes. As HIV specifically acts on the immune system, 
destroying or impairing its function, the most common causes 
of death for HIV patients are various pathogenic infections 
(including fungal, viral, bacterial, and parasitic infections) and 
cancer (Bonnet et al., 2005; Antiretroviral Therapy Cohort 
Collaboration, 2010; Palladino et al., 2011). 

Both parasites and HIV infection have major effects on 
the host immune system, and co-infection is widespread 
(Bentwich, 2000). Helminths can induce a Th2 bias, leading to 
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suppression of Th1 responses to HIV, as well as to expansion 
of Th2 lymphocytes that are more susceptible to HIV infection 
(Bentwich et al., 1999). Several epidemiological and systematic 
studies claim that deworming is a preventive and possibly 
therapeutic measure in patients co-infected with helminths 
and HIV (Borkow & Bentwich, 2006; Walson et al., 2009; 
Walson & John-Stewart, 2007), thus indicating the potential 
side effect of parasites. However, a more recent meta-analysis 
showed that deworming has little or no favorable effect on 
HIV treatment (Means et al., 2016). Unfortunately, direct 
experimental evidence from animal models is still missing, and 
the relationship between parasites and HIV remains complex 
and unclear. 

Non-human primates (NHPs) are crucial animal models for 
understanding the pathogenesis of HIV infection (Karlsson 
et al., 1997; Letvin, 1992). Simian immunodeficiency virus 
(SIV) infections in rhesus macaques and northern pig-tailed 
macaques (NPMs) result in an immunocompromised and 
AIDS-like disease similar to that seen in HIV patients, as 
established in our laboratory (Pang et al., 2017; Zhang et 
al. 2017a, b). In this study, NPMs with different parasitic 
loads after infection with SlIVmac239 were used to detect 
immunological and pathological differences to reveal the 
relationship between viral infection and parasites. 


MATERIALS AND METHODS 


Ethics statement 

All animal experiments were performed according to the 
guidelines of the Committee on Animals of the Kunming 
Institute of Zoology, Chinese Academy of Sciences (approval 
number: SYDW-2015023). NPMs included in this study were 
housed at the Kunming Primate Research Center under the 
guidance of the American Association for Assessment and 
Accreditation of Laboratory Animal Care (AAALAC). 


Animals and sample collection 

After evaluating each individual for B virus, SIV, STLV, and SRV 
by PCR, six negative male NPMs (5-7 years old, numbered 
08247, 08287, 09203, 09211, 10205, and 10225) were enrolled 
in this study. Each NPM was inoculated with 2 mL of 
normal saline containing a 3 000 50% tissue culture infection 


dose (TCIDso) of SlVmac239 into the posterior tibial vein. 


Peripheral blood was collected into an EDTAK» vacutainer by 
venipuncture. Subsequently, peripheral blood mononuclear 
cells (PBMCs) were separated by Ficoll density gradient 
centrifugation (2 500 r/min, 20 min). PBMCs and plasma were 


stored in liquid nitrogen and frozen at —80 ?C , respectively. 


Ketamine hydrochloride (0.2 mL/kg) was used to anaesthetize 
animals prior to experiments. 


Parasite detection 

Pre-infection and at the 50th week after viral infection, plasma 
from the NPMs was collected and evaluated for parasitic 
status. In total, 14 types of parasites were included in 
this study. Entamoeba histolytica was detected by thick 


smear, as described previously (Warhurst & Williams, 1996). 


Trypanosoma brucei, Cryptosporidium parvum, Plasmodium 


knowlesi, Toxoplasma gondii, Giardia lamblia, Echinococcus 
granulosus, Clonorchis sinensis, Schistosoma japonicum, 
Paragonimus spp., filarial worms, and Trichinella spiralis were 
detected by their respective parasitic antigen ELISA kits 
(MEIMIAN, China) following the manufacturer's instructions. 
Spirometra mansoni and Angiostrongylus cantonensis were 
detected by their respective human antibodies ELISA 
kits (MEIMIAN, China; JIANLUN, China) following the 
manufacturer’s instructions. For all ELISA experiments, the 
Calculate Critical (CUT OFF) value was equal to the average of 
the negative control well plus 0.15. The well with an OD greater 
than or equal to the CUT OFF was regarded as a positive result. 


Full blood differential counts 

Before separation of PBMCs, 80 uL of blood was removed from 
the EDTAKs vacutainer. Whole blood differential counts were 
performed by an auto hematology analyzer (RR-OB101338, 
MINDRAY, Shenzhen, China) according to the manufacturer's 
instructions. 


Plasma viral loads 

Viral loads in plasma were measured by quantitative 
PCR for SIV gag RNA, as described previously (Tian et 
al, 2015) (SIV F: 5" TCGGTCTTAGCTCCATTAGTGCO-3'; 
SIV R: 5'-GCTTCCTCAGTGTGTTTCACTTTC-3'; SIV probe: 
5'-CTTCTGCGTGAATGCACCAGATGACGC-3’). 


Flow cytometry 

Multicolor flow cytometric analysis was performed as described 
previously (Xia et al., 2009; Zheng et al., 2014). Anti-CD3 PE 
(clone SP34-2), anti-CD8 PECy7 (clone RPA-T8), anti-CD20 
FITC (clone 2H7), anti-CD14 APC (clone M5E2), anti-CD3 
APC-Cy7 (clone SP34-2), and anti-HLA-DR APC (clone G46-6) 
were purchased from BD Pharmingen (Franklin Lakes, New 
Jersey, USA). Anti-CD4 PerCP (clone OKT4) was purchased 
from Biolegend (San Diego, CA, USA). Anti-CD38 FITC 
(clone AT-1) was purchased from STEMCELL (Vancouver, 
Canada). Anti-PD-1 PE (clone eBioJ105) was purchased from 
eBioscience (San Diego, CA, USA). Samples were processed 
on a BD FACSVerse flow cytometer (Franklin Lakes, New 
Jersey, USA), and data were analysed using FlowJo software 
(vX.0.7, Tree Star). 


Tissue collection and histopathology 

Both 08247 and 09211 died naturally on the 70th and 68th 
weeks after viral infection, respectively. Tissues were fixed 
in 1096 formaldehyde and later embedded in paraffin. The 
paraffin-embedded samples were cut into 4-um-thick sections 
and stained with hematoxylin-eosin (H&E). All samples were 
photographed and examined using a LEICA DMI 4000B 
microscope (Germany). 


Statistical analysis 

For all data, GraphPad 6.0 was used to perform statistical 
analyses. Two-way analysis of variance (ANOVA) was used 
to compare all data. P<0.05 were regarded as statistically 
significant. 
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RESULTS 


Parasitic infection status and categories 

The parasitic infection statuses of the NPMs are shown in 
Table 1. None of the NPMs was infected with Plasmodium 
knowlesi, Echinococcus granulosus, Clonorchis sinensis, 
Paragonimus spp., Angiostrongylus cantonensis, Trichinella 
spiralis, or Spirometra mansoni. Both 08247 and 08287 
showed positive parasitic reactions before viral infection, and 
negative reactions after infection. However, 09203, 09211, 
10205, and 10225 exhibited substantially different results. At 
the 50th week after infection, the parasites detected before viral 
infection remained, and no new parasites were detected after 
viral infection. Therefore, we grouped 08247 and 08287 into 
the Parasites Exit (PE) group and the other four NPMs into 
the Parasites Remain (PR) group to evaluate the relationship 
between SIV and parasitic infection. 


Higher viral loads in the PE group 

NPMs in both groups had peak plasma viral loads at the second 
to third week post-infection, after which the viral loads showed 
a downward trend until the sixth week after injection (Figure 
1). Afterwards, the PE and PR groups showed substantially 
different trends. Viral loads in the PE group stopped 


decreasing and generally stabilized at 10*—105 copies/mL. 


However, the PR group showed a continuing downward trend, 
and approached detection limits at the 40th—50th weeks after 
infection. In this study, the 50 experimental weeks were divided 
into SIV acute phase and SIV chronic phase at the 12th week 
post infection to perform additional analysis (Pantaleo et al., 
1993b). Compared with the PR group, the PE group showed 
significantly higher plasma viral loads in the chronic stage 
(P<0.000 1). 


Lower CD4* T cell levels in the PE group 

Absolute counts of CD4* T cells, the primary cells infected by 
HIV or SIV, were evaluated by flow cytometry. Compared with 
the PR group, NPMs in the PE group had significantly fewer 
CD4* T cells in peripheral blood (P«0.000 1, Figure 2A). This 
was in accordance with the results in Figure 2B, which show 
a significantly lower CD4*/CD8* T cell rate in the PE group 
(P«0.001). These two results strongly indicated that the NPMs 
in the PE group had more seriously damaged CD4* T cells. 


Immune activation and exhaustion 

Chronic immune activation is proposed to be a key determinant 
of AIDS pathogenesis (Rajasuriar et al., 2013), and was 
marked by CD38 and HLA-DR in this study. There were 
significant differences between the PE and PR groups in CD38 
expression during the whole study period and in HLA-DR 
expression during the chronic phase in CD4* and CD8* T cells 
(Figure 3A, B, D, E). These data strongly indicated that NPMs 
in the PE group showed more obvious immune activation than 
those in the PR group during SIV infection. Additionally, the PE 
group had a higher level of PD-1 expression in CD4* T cells 
during the whole study period (P«0.000 1) and in CD8* T cells 
during the chronic phase (P=0.027 3), as shown in Figure 3C 
and 3F respectively. Collectively, these results revealed that 
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NPMs in the PE group showed higher immune activation and 
more serious exhaustion of CD4* and CD8* T cells. 


Cell distribution in peripheral blood 

The fluctuant cell distribution in peripheral blood was shown in 
Figure 4. The number of white blood cells showed no significant 
differences between the two groups (Figure 4C). However, in 
the SIV acute and chronic phases, NPMs in the PE group 
showed fewer red blood cells compared with the PR group, as 
shown in Figure 4A. The number of platelets also significantly 
declined in the PE group compared with the PR group during 
the SIV chronic phase (Figure 4B), which is important given 
that platelet count is a main marker for immunocompromised 
thrombocytopenic purpura induced by HIV or SIV (O'Bryan et 
al., 2015). 


Dynamics of monocytes, B cells, and T cells in peripheral 
blood 

Cell counts before viral infection were calculated as 10096. The 
dynamics of the three cell types were evaluated, as shown 
in Figure 5. Monocytes, which play an important role in 
innate immunity, showed a significant increase in the PE group 
compared with the PR group (Figure 5A). Conversely, T cells 
and B cells, markers of acquired immunity, were significantly 
reduced in the PE group compared with the PR group, as 
shown in Figure 5B, C. 


Pathological observation of injured organs 

In total, 22 organs from each NPM were collected and stained 
by H&E, four of which displayed obvious damage, as shown 
in Figure 6. Compared to PR-09211, which maintained almost 
normal hepatic lobule appearance, liver damage in PE-08247 
was obvious (Figure 6A). Alimentary canal damage is strongly 
related to chronic immune activation (Griffin, 1990; Mestecky 
et al., 2009), and was detected in our study. In Figure 6C, D, 
obvious damage in the ileum of both 08247 and 09211 could 
be seen. Compared to 09211, PE-08247 demonstrated more 
serious injuries in the cecum without any observed enteraden 
and in the colon with a discrete and defective mucosal surface. 


Pathological observation of immune organs 

The spleen has no direct connection to the lymphatic system; 
instead, it collects antigens from the blood and is involved 
in immune responses to blood-borne pathogens (Murphy & 
Weaver, 2016). Compared to 09211, the white pulp of 
08247 in the PE group had severe fibrosis and decreased 
lymphocytes, as shown in Figure 7A. Lymph nodes are the 
most important organ in HIV or SIV infection, spread, and 
damage to lymphocytes (Pantaleo et al., 1993a). As shown in 
Figure 7C—F, lymphoid follicles could be observed in all images, 
and there were no significant differences between the PE and 
PR groups in axillary and inguinal lymph nodes. However, 
compared to 09211 in the PR group, mesenteric lymph nodes 
in 08247 in the PE group were observed with lymphoid follicle 
fibrosis and decreased lymphocytes in the paracortex, strongly 
indicating more serious damage in the PE group (Figure 7G). 


Table 1 Parasitic infection status 
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Figure 2 Dynamics of CD4* T cells in peripheral blood 
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Figure 3 Expression of immune activation and exhaustion markers 
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Figure 4 Cell distribution in peripheral blood 
A: Red blood cell (RBC) counts; B: Platelet (PLT) counts; C: White blood cell (WBC) counts. 
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Figure 5 Dynamics of monocytes, B cells, and T cells in peripheral blood 
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Figure 6 Histological changes in liver, ileum, cecum, and colon 
A, B: Liver (x200); C, D: Ileum (x100); E, F: Cecum (x 100); G, H: Colon (x 100). Red arrow: Fibrosis; Blue arrow: Hemorrhage; Yellow arrow: Hepatic lobule; 
Black arrow: Enteraden; White arrow: Connective tissue exposure; Red rectangle: Mucosal surface. 
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Figure 7 Histological changes in immune organs (x 100) 
A, B: Spleen; C, D: Axillary lymph nodes; E, F: Inguinal lymph nodes; G, H: Mesenteric lymph nodes. Red arrow: Abnormal white pulp; Yellow arrow: White pulp; 


Black arrow: Lymphoid follicle; White arrow: Abnormal lymphoid follicle. 
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DISCUSSION 


In this study, no parasites were detected in 08247 and 08287 
in the PE group at the 50th week after viral infection, although 
both were positive for parasites before viral injection. — In 
addition, in the PE group, more significant damage to the 
mucosal barrier, mesenteric lymph nodes, and immune cells 
were observed. 


CD4* T cell counts and viral loads are the most important 
indices used to predict disease progress and immune system 
damage (Rowley, 2014). The persistent decrease in CD4* 
T cell counts with viral load increase represents serious 
immunological compromise and poor prognosis, which was 
observed in the PE group. In addition, immune activation 
is a sensitive index for predicting disease progression and is 
central to the understanding of HIV-1 disease pathogenesis 
and progression in association with inflammation (lpp et al., 
2014). Therefore, CD38 and HLA-DR, two markers of immune 
activation, were detected by flow cytometry in this study, which 
indicated more significant immune activation in the PE group 
compared with the PR group. 


The loss of enteric mucosal barrier integrity is associated 
with rapid depletion of enteric lymphocytes, immune activation, 
and microbial translocation (Brenchley et al., 2006; Douek, 
2007). As shown in Figure 6, PE-08247 exhibited serious 
damage to the mucosal barrier. Unsurprisingly, antigen origin 
from the translocation of microbes and other bacterial products 
stimulates mesenteric lymph nodes, resulting in lymphoid 
tissue hyperplasia and ultimately in the diffuse effacement of 
lymph node architecture (Lederman & Margolis, 2008), which 
was observed in the PE group (Figure 7). The decrease in 
platelets, which is correlated with immunological compromise 
and virus-related mechanisms, was also detected in the PE 
group. The significant increase in monocytes and decrease 
in T cells in peripheral blood indicated changes to the immune 
system in the PE group. These results showed that NPMs in 
the PE group were more seriously immunocompromised. 


Considering the parasitic detection results in the PE 
group before and after viral infection, parasites may exit 
hosts that are seriously immunocompromised. Doenhoff et 
al. (1986) found that granuloma formation and consequent 
schistosome egg excretion might be reduced in animal models 
of immunosuppression. This was supported by subsequent 
studies in humans (Karanja et al., 1997; Mwanakasale et 
al., 2003; N’Zoukoudi-N’Doundou et al., 1995). Conversely, 
another hypothesis suggests that schistosomal infection could 
be regarded as an immune reconstitution marker following 
antiretroviral therapy, whereby pathological responses to 
schistosomes and schistosome eggs are inhibited in advanced 
HIV disease and deteriorate with immune recovery (De Silva et 
al., 2006; Fernando & Miller, 2002). Additionally, in contrast 
to the theory that deworming is a preventive and possible 
therapeutic measure for patients with co-infection of helminths 
and HIV (Borkow & Bentwich, 2006; Walson et al., 2009; 
Walson & John-Stewart, 2007), recent meta-analysis provides 
stronger evidence that deworming has little or no favorable 


effect on HIV treatment (Means et al., 2016). These studies 
suggest that changes in the host induced by HIV or SIV 
produce adverse effects on parasites, and parasites might not 
contribute to the progression of HIV disease, matching the 
hypothesis of Brown et al. (2006). 

There were three notable limitations to this study. First, there 
were only six NPMs enrolled, which might not fully represent 
the differences between PE and PR groups. Second, feces 
detection is regarded as the gold standard for the detection of 
some parasites, particularly intestinal parasites. Unfortunately, 
we did not collect fecal samples from the six NPMs before 
viral infection. Therefore, all experiments were performed 
using antigen or antibody ELISA kits. Finally, tissues were 
collected from only two NPMs. The other survivors were 
kept for long-term parasitic experiments. Euthanasia was not 
performed in this study. However, most of the results in this 
study showed an identical trend suggesting the conclusions 
were effective and meaningful. 


CONCLUSIONS 


To the best of our knowledge, this study is the first to provide 
direct evidence on the changes in parasitic infection status 
before and after SIV injection in an animal model. Both HIV 
and SIV destroy the host immune system, leading to serious 
immunological compromise, which enhances the possibility 
for pathogenic microorganisms, including bacteria, fungi, and 
viruses, to invade and cause damage. In contrast, our results 
indicated that parasites may exit an immunocompromised host. 
These prospective results should be evaluated further with a 
larger sample size. 
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ABSTRACT 


Globally, coxsackievirus B4 (CV-B4) has been 
continuously isolated and evidence suggests an 
association with the development of pancreatitis and 
type | diabetes. In addition, CV-B4 is also associated 
with myocarditis and severe central nervous system 
(CNS) complications, which remain poorly studied 
and understood. In the present study, we established 
an Institute for Cancer Research (ICR) mouse model 
of CV-B4 infection and examined whether CV-B4 
infection resulted in a predisposition to myocarditis 
and CNS infection. We found high survival in both 
the treatment and control group, with no significant 
differences in clinical outcomes observed. However, 
pathological lesions were evident in both brain and 
heart tissue of the CV-B4-infected mice. In addition, 
high viral loads were found in the neural and cardiac 
tissues as early as 2 days post infection. Expressions 
of IFN-y and IL-6 in sera were significantly higher in 
CV-B4-infected mice compared to uninfected 
negative controls, suggesting the involvement of 
these cytokines in the development of histopathological 
lesions. Our murine model successfully reproduced 
the acute myocarditis and cerebral cortical neuron 
edema induced by CV-B4, and may be useful for the 
evaluation of vaccine candidates and potential 
antivirals against CV-B4 infection. 


Keywords: Coxsackievirus B4; Myocarditis; CNS; 
Edema; Neurons 


INTRODUCTION 


Coxsackievirus B (CV-B) comprises a number of established 
cardiomyopathy-associated viral serotypes, taxonomically 
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classified within the genus Enterovirus species B, family 
Picornaviridae. There are currently 63 recognized enterovirus 
serotypes (www.picornaviridae.com) with diverse tissue tropisms 
and correspondingly broad spectrum of disease presentation, 
including myocarditis, encephalitis, paralysis, meningitis, upper 
and lower respiratory disease, pleurodynia, herpangina, 
myopericarditis, pancreatitis, and type | diabetes (Knowles et 
al., 2011; Tapparel et al., 2013). In the prior two decades, the 
relationship between CV-B and myocarditis/dilated cardiomyopathy 
has been substantiated. For example, CV-B infections have 
been identified in approximately 2596—4096 cases of acute 
myocarditis and dilated cardiomyopathy in young adolescents 
and adults (Gaaloul et al., 2014; Leonard, 2004). Furthermore, 
an epidemiological survey has revealed the co-circulation of six 
serotypes of CV-B in a single location in Changchun, Jilin 
Province, China, with the predominant agents of viral 
myocarditis in infants from 2002-2011 being CV-B3 (30.196— 
36.5%), CV-B4 (18.096—24.896), and CV-B5 (3.396—8.196) 
(Zhang et al., 2014). Significantly, CV-B4 is also associated with 
severe central nervous system (CNS) complications and 
neurological sequelae (Zhu et al., 2015). 

Previous research has employed mouse models to study CV- 
B4 infection. These studies have mainly focused on the 
inflammatory insults in exocrine tissue damage (De Palma et 
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al., 2009) and the association with type 1 diabetes (McCall et 
al., 2015) during CV-B4 infection. Due to the lack of a stable 
mouse model of CV-B4 infection, the exact mechanism(s) of 
how CV-B4 results in myocarditis and CNS invasion remain 
unclear. Furthermore, pathogenesis studies are impeded by the 
absence of appropriate animal models of disease that would 
facilitate the development of prophylactic vaccines and antivirals 
for the prevention and treatment of acute disease. 

We investigated whether CV-B4 infection results in a 
predisposition to myocarditis and CNS infection in Institute for 
Cancer Research (ICR) neonatal mice. Dynamic detection of viral 
loads, hematoxylin and eosin staining, and immunohistochemical 
examination were performed to determine tissue tropisms, 
pathological lesions, and distribution of CV-B4 in vivo, respectively. 
Finally, the expression levels of virus-induced inflammatory cytokines 
were examined and correlated with viral titers and clinical scores. 


MATERIALS AND METHODS 


Virus and cells 

Human laryngeal carcinoma epithelium (Hep-2) cells were 
cultured in Dulbecco's modified Eagle's medium (DMEM; 
Gibco, USA), supplemented with 10% fetal bovine serum 
(Gibco) and 1% penicillin-streptomycin at 37 °C under a 5% 
CO» atmosphere. The CV-B4 LY114F strain employed in the 
present study was isolated from the stool of a three-year-old 
patient who presented with classical symptoms of hand, foot, 
and mouth disease (HFMD) in Linyi city, Shandong Province, 
China, in 2015. Virus propagation in Hep-2 cells and 5096 tissue 
culture infective doses (TCIDso) were determined in accordance 
with the methods of Reed & Muench (1938). 


Ethics statement and animal infection model 

Institute for Cancer Research (ICR) mice were purchased from 
Beijing Vital River Laboratory Animal Technology Co., Ltd. The 
animal experiments were approved by the Taishan Medical 
College Administrative Committee for Laboratory Animals, and 
all procedures involving animals were performed in accordance 
with the Shandong Laboratory Animal Welfare and Ethics 
Administrative Committee. 

Three-day-old ICR mice were challenged with the CV-B4 
LY114F strain (107 TCIDayanimal) by intramuscular (i.m.) injection, 
with control groups injected with DMEM. Body weights, clinical 
manifestations, and survival rates of mice infected with LY114F 
were monitored daily until 10 days post infection (dpi). Multivariate 
analysis was used to analyze the differences in body weights 
between infected mice and negative controls at each observation 
time point. Clinical criteria were scored as follows: 0, healthy; 1, 
lethargy and inactivity; 2, hind limb weakness; 3, single limb 
paralysis; 4, double hind limb paralysis; and 5, death or dying. 


Dynamics of CV-B4 viral RNA titers in infected mice 

Brain, heart, contralateral hind limb skeletal muscle, lung, 
intestine, spleen, and blood samples from each mouse (n=3 per 
time point) infected with LY114F were collected at 1—5 dpi, 
respectively. Viral RNA was extracted from equivalent weights 
of tissue samples and volumes of blood samples from infected 


and control mice with TRIzol reagent (Takara, China), as per 
the manufacturer's instructions. cDNA was generated using a 
reverse transcription kit (Takara, China) for 45 min at 42 °C and 
a GoldStar TaqMan Mixture kit (CWBIO) was employed for real- 
time PCR. Oligonucleotide primers and hydrolysis probe 
sequences are as follows: sense, 5-CCTGAATGCGGCTAAT 
CC-3'; antisense, 5'-TTGTCACCATWAGCAGYCA-3’; and hydrolysis 
probe, 5'-FAM-CCGACTACTTT GGGWGTCCGTGT-BHQ1-3'. 
Real-time PCR thermocycling was performed on a LightCycler 
96 platform with the following conditions: 95 °C for 10 min; 40 
cycles of 95 °C for 15 s, and 60 °C for 1 min. Standard curves 
for absolute quantification of CV-B4 copy numbers in different 
tissue and blood samples were established, as described 
previously (Zhang et al., 20172, b). 


ELISA measurement of cytokine levels in sera 

Peripheral blood was collected from 3-day-old mock-infected 
controls or mice infected with a 10’ TClDso/animal dose of 
LY114F from 1—5 dpi (n=3 per time point). Peripheral blood was 
centrifuged immediately at 10 000 r/min for 10 min at room 
temperature and sera were stored at —80 ?C prior to use. The 
levels of IFN-y, IL-6, IL-4, IL-18, IL-10, IL-13, IL-18, and TNF-a in 
sera were determined by commercial ELISA kits (Multisciences 
Biotechnology, China), according to the manufacturer's protocols. 


Histopathological and immunohistochemical staining 
Three-day-old mice were infected with LY114F (107 TCIDso/ 
animal). Brains, hearts, muscles, and lungs from experimental 
mice and controls were subjected to histopathological and 
immunohistochemical examination. Tissues were fixed in 1096 
neutral buffered formalin for 72 h and then embedded in 
paraffin, as per previously described procedures (Yu et al., 
2014). For pathological examinations, formalin-fixed, paraffin- 
embedded sections (4 um thick) stained with hematoxylin and 
eosin were employed. For immunohistochemical examination, 
paraffin-embedded sections were dewaxed, dehydrated, and 
microwaved for 15 min at 95 °C to 99 °C in EDTA buffer (1 
mmol/L). Polyclonal mouse anti-CV-B4 antibody (1:300 dilution; 
Abcam, UK) was applied for 15 h at 4 °C. A secondary 
horseradish peroxidase-conjugated goat anti-mouse IgG (1:400 
dilution; ZSGB-BIO, China) was then applied for 30 min at room 
temperature, followed by avidin-biotin-peroxidase complex and 
3,3'-diaminobenzidine tetrahydrochloride chromogen (1:1000 
dilution; ZSGB-BIO, China). Subsequently, tissue sections were 
counterstained with hematoxylin and negative controls were 
incubated with PBS instead of primary antibody. 


Statistical analysis 
All statistical analyses were performed with the Statistical Analysis 
System 9.2 (SAS, USA). Differences in mean tissue viral load 
and serum cytokine concentrations were determined with two- 
tailed analysis of variance (ANOVA). A difference was considered 
significant at P«0.05. 


RESULTS 


Establishment of the CV-B4 infection mouse model 
Multivariate analysis showed that the body weights of the 
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infected mice were significantly higher than those of the 
negative controls at 6 dpi (F=5.439, P=0.031), 8 dpi (F=9.521, 
P=0.008), 9 dpi (F=8.526, P=0.009), and 10 dpi (F=10.601, 
P=0.004) (Figure 1). There were no significant differences 
between the CV-B4 LY114F-treated mice and the control 
group with respect to mean clinical score, and survival rate. 
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Figure 1 Body weights of ICR mice infected with the CV-B4 LY114F 
strain 

Control animals were inoculated with saline. All mice were monitored 
daily for body weight until 10 dpi. 


Histopathological and immunohistochemical examination 
To examine the histopathological lesions and distribution of viral 
antigens in mice post infection, histopathological analysis and 
immunohistochemical examination of the brains, hearts, 
Skeletal muscles, and lungs were performed. Results showed 
evident pathological lesions in the brain (Figure 2A-D) and 
heart (Figure 2E-H) tissues in mice infected with CV-B4 
LY114F, including cerebral cortical neuron edema (Figure 2B) 
and myocardial lymphocytic infiltration (Figure — 2F). 
Furthermore, immunohistochemical analysis demonstrated the 
presence of CV-B4 antigens in brain neurons (Figure 2D) and 
myocardial fibers (Figure 2H); however, there were no distinct 
histopathological lesions or antigen detectable in the lung 
(Figure 2I-L) or skeletal muscle tissues (Figure 2M—P). 


Dynamics of CV-B4 viral RNA titers in the CV-B4 murine 
model 

Viral loads in 3-day-old mice inoculated with 10" TCIDso of CV- 
B4 LY114F were measured at 1—5 dpi. The kinetics of viral 
loads in the brain, heart, contralateral hind limb skeletal muscle, 
lung, intestine, spleen, and blood samples are illustrated in 
Figure 3. In the infected group, CV-B4 RNA was detectable in 
the brain at 3 dpi and peaked at 5 dpi with a viral load of >10° 
copies/mg tissue. Notably, viral loads in the heart were the 
earliest evident of all organs tested (2 dpi) and gradually 
decreased from 2.4x10’ copies/mg (2 dpi) to 1.1x10" copies/ 
mg (5 dpi). Additionally, there was a slow increase in the viral 
titer following infection from days 1 to 5 in lung tissues and 
skeletal muscle, in which viral loads each peaked at 
approximately 3.0x10" copies/mg at 5 dpi. In contrast, from 1 to 
5 dpi, viral titers in the blood (0.0—1.0x109 copies/mL), intestine 
(0.0—3.1x109 copies/mg), and spleen (0.0—3.3x109 copies/mg) 
were much lower. 
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Expression of inflammatory cytokines in sera 

The expressions of inflammatory cytokines, including IFN-y, IL- 
6, IL-10, IL-13, TNF-a, IL-18, IL-48, and IL-4, in sera were 
measured by ELISA at 1 to 5 dpi in 3-day-old ICR mice infected 
with CV-B4 LY114F (10’ TClIDao/animal dose). The expressions 
of IFN-y and IL-6 in sera were significantly higher in CV-B4- 
infected mice compared to uninfected negative controls. In CV- 
B4-infected animals, the expression level of IFN-y (Figure 4A) 
in serum peaked very early (~7 100 pg/mL at 1 dpi), and then 
decreased rapidly to «4 000 pg/mL by 2 dpi; however, high 
expression levels (>2 000 pg/mL) were maintained up to 5 dpi. 
In contrast, the expression level of IL-6 (Figure 4B) in serum 
was low at 1 dpi, but increased steadily during infection and 
peaked at 5 dpi with a concentration of >500 pg/mL. The titers 
of IL-10 (Figure 4C) were detectable early in infection (10 pg/mL 
at 1 and 2 dpi) compared to uninfected controls and peaked at 
16.44 pg/mL by 3 dpi, with low expression levels maintained 
throughout the infection. The titers of IL-13 (Figure 4D) were 
also detectable early in infection and peaked at 20.11 pg/mL at 
2 dpi, before gradually decreasing. The expression of TNF-a 
(Figure 4E) was lower at earlier time points, and peaked at 
17.37 pg/mL by 4 dpi. The titers of IL-18 (Figure 4F) in both 
experimental and control groups showed high levels of 
expression (72100 pg/mL) across all time points, and there were 
no significant differences between the groups. Finally, the 
expressions of IL-18 and IL-4 were not detectable in either 
group. 


DISCUSSION 


Previous studies on the CV-B4 viral serotype have primarily 
focused on cardiotropism, with little emphasis on the 
pathogenicity of CV-B4 in the CNS (El Hiar et al., 2012; Hu et 
al., 2012). Earlier reports demonstrated that >90% of reported 
CV-B infections occur in infants, children, and adolescents 
(Romero, 2008). Therefore, establishing a sensitive neonatal 
animal model of CV-B4 infection is important for pathogenesis 
studies of CV-B4-induced CNS disease in infants. Such models 
could clarify the pathological processes and enable 
development of vaccines and antivirals for the prevention and 
treatment of acute neurological cases. 

After 3-day-old neonatal mice were inoculated with 10" 
TCIDso CV-B4 in the present study, viral RNA was detectable in 
several organs, including brain, heart, lung, and skeletal 
muscle. Of note, viral titers in the brain were three-fold greater 
than that in the other organs at the latter stage of infection (5 
dpi). In addition, readily discernible pathological changes were 
detected in both brain and myocardial tissue. Compared with 
enterovirus species A virus infection models, including EV-A71 
(Dong et al., 2016; Yue et al., 2016), CV-A16 (Huang et al., 
2015; Mao et al., 2012), CV-A6 (Zhang et al., 2017b), and CV- 
A10 (Zhang et al., 20178), which are associated with hind limb 
paralysis and/or severe respiratory involvement, the enterovirus 
species B CV-B4 had comparably minor pathology in skeletal 
muscle and lung tissue. 





Figure 2 HE and IHC analyses of infected 3-day-old mice after i.m. challenge with the CV-B4 LY114F strain 

No histological changes were observed in the brain, heart, lung, or hind limb muscle tissue of the negative control mice (A, E, I, and M), as well as in the 
lung (J) or muscle (N) tissues of the infected mice. Infected mice exhibited cerebral cortical neuron edema in the brain (B) and myocardial lymphocytic 
infiltration in heart tissue sections (F), as indicated by arrows. IHC analyses showed the existence of viral antigen in the brain and heart tissue of 
infected mice (D, H). There was no antigen detected in the lung (L) and muscle (P) tissue of the infected mice and in the mice challenged with the PBS 
buffer (C, G, K, and O). Magnification: x200, scale bars: 100 um. 
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Figure 3 Mean viral loads in tissue and peripheral blood of CV-B4-infected mice 

Three-day-old ICR mice were i.m. inoculated with 107 TCIDso of the CV-B4 LY114F strain. Viral loads (109 copies/mg or mL) in the brain, heart, blood, 
lung, intestine, spleen, and muscle of infected mice were quantified by RT-qPCR. Data represent the mean results of three mice + standard deviation of 
the mean. 
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Figure 4 Peripheral blood cytokine expression levels in neonatal mice infected with the CV-B4 LY114F strain 

Levels of IFN-y (A), IL-6 (B), IL-10 (C), IL-13 (D), TNF-a (E), and IL-18 (F) in the sera of 3-day-old ICR mice i.m. inoculated with LY114F and of negative 
controls (NC) at 1, 2, 3, 4, and 5 dpi were determined using murine ELISA assays. Expression levels of cytokines in sera between experimental and 
control groups were compared by Fisher's least significant difference multiple comparison test. *: P<0.01. Data represent the mean results of three mice+ 


standard deviation of the mean. 


The interaction between pathogen and the host immune system 
is critical to develop an understanding of the mechanisms of 
viral pathogenesis. Clinical studies have shown that the high 
expressions of several cytokines, such as IL-6, TNF-a, and IL- 
10, are closely correlated with clinical manifestations and 
complications associated with CV-B4 infection (Alidjinou et al., 
2013; Gu et al., 2009). Our data showed that after neonatal 
mice were inoculated with CV-B4, there was a gradual increase 
in the level of IL-6 in serum, different from that observed in 
enterovirus species A virus infection (CV-A6 and CV-A10) in 
ICR neonatal mice (Zhang et al., 2017a, b). Our previous 
studies on CV-A6 and CV-A10 mouse models of infection also 
showed high expression levels of IL-6 throughout all stages of 
infection (Zhang et al., 2017a, b). Therefore, we speculate that 
the differences in pathogenicity between CV-A and CV-B 
enteroviral species in the same neonatal ICR murine model 
may be attributable to differences in skeletal muscle tropisms 
and proinflammatory cytokine expression levels, such as that of 
IL-6, in peripheral blood. Previous in vivo experiments studying 
EV-A71 infection have shown that treatment with anti-IL-6 
neutralizing antibodies increases survival rates, whereas 
treatment with IL-6 can exacerbate pulmonary dysfunction in 
EV-A71-infected mice (Khong et al., 2011). We also observed 
this latter phenomenon of disease exacerbation in CV-A6 and 
CV-A10 murine models of infection (Zhang et al., 2017a, b). 
The coxsackievirus and adenovirus receptor (CAR) is 
responsible for coxsackie B virus infection in human cells 
(Bergelson et al., 1997), with non-permissive cells transfected 
with mCAR cDNA rendered susceptible to CV-B viral infections 
(Bergelson et al., 1997; Bergelson et al., 1998). CAR is reported 
to be employed by laboratory reference strains and clinical 
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isolates of all six serotypes of EV-B (Martino et al., 2000); 
however, different clinical CV-B isolates have been found to 
possess distinct interactions with CAR (Riabi et al., 2014). In 
addition, enhanced CAR expression is associated with 
experimental autoimmune myocarditis in adult mice (Ito et al., 
2000) and treatment of CV-B3-infected BALB/c mice with 
CAR4/7 can aggravate cardiac injury (Dörner et al., 2006). 
Although CAR expression is significantly higher in dilated 
cardiomyopathy (DCM) cases than in negative controls, no 
significant differences in EV viral loads between DCM and non- 
DCM cases have been observed (Sharma et al., 2016). 
Consequently, the relationship between increased CAR 
expression and viral load of enteroviruses, including CV-A and 
CV-B species, remains unclear and warrants further investigation. 

In conclusion, ICR neonatal mice i.m. infected with CV-B4 
exhibited significant pathology in the brain and myocardium. 
This represents a suitable model for the establishment of CV-B4 
infection to study the pathogenesis of CNS and cardiac 
complications and is a resource for the development of both 
prophylactic vaccines and antivirals to reduce morbidity and 
mortality in children. 
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